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We study the role of particle-hole symmetry on the universality class of various quantum phase 
transitions corresponding to the onset of superfluidity at zero temperature of bosons in a quenched 
random medium. To obtain a model with an exact particle-hole symmetry it is necessary to use 
the Josephson junction array, or quantum rotor, Hamiltonian, which may include disorder in both 
the site energies and the Josephson couplings between wavefunction phase operators at different 
sites. The functional integral formulation of this problem in d spatial dimensions yields a (d -I- 1)- 
dimensional classical XK-model with extended disorder, constant along the extra imaginary time 
dimension — the so-called random rod problem. Particle-hole symmetry may then be broken by 
adding nonzero site energies, which may be uniform or site-dependent. We may distinguish three 
cases: (i) exact particle-hole symmetry, in which the site energies all vanish, (ii) statistical particle- 
hole symmetry in which the site energy distribution is symmetric about zero, vanishing on average, 
and (iii) complete absence of particle-hole symmetry in which the distribution is generic. We explore 
in each case the nature of the excitations in the non-superfluid Mott insulating and Bose glass phases. 
We show, in particular, that, since the boundary of the Mott phase can be derived exactly in terms 
of that for the pure, non-disordered system, that there can be no direct Mott-superfluid transition. 
Recent Monte Carlo data to the contrary can be explained in terms of rare region effects that 
are inaccessible to finite systems. We find also that the Bose glass compressibility, which has the 
interpretation of a temporal spin stiffness or superfluid density, is positive in cases (ii) and (iii), but 
that it vanishes with an essential singularity as full particle-hole symmetry is restored. We then 
focus on the critical point and discuss the relevance of type (ii) particle-hole symmetry breaking 
perturbations to the random rod critical behavior, identifying a nontrivial crossover exponent. This 
exponent cannot be calculated exactly but is argued to be positive and the perturbation therefore 
relevant. We argue next that a perturbation of type (iii) is irrelevant to the resulting type (ii) critical 
behavior: the statistical symmetry is restored on large scales close to the critical point, and case (ii) 
therefore describes the dirty boson fixed point. Using various duality transformations we verify all 
of these ideas in one dimension. To study higher dimensions we attempt, with partial success, to 
generalize the Dorogovtsev-Cardy-Boyanovsky double epsilon expansion technique to this problem. 
We find that when the dimension of time < e!; ~ ^ is sufficiently small a type (ii) symmetry 
breaking perturbation is irrelevant, but that for sufficiently large e-r > e% particle-hole asymmetry 
is a relevant perturbation and a new stable fixed point appears. Furthermore, for Et > e^'^ ~ |, this 
fixed point is stable also to perturbations of type (iii): at e = e^^ the generic type (iii) fixed point 
merges with the new fixed point. We speculate, therefore, that this new fixed point becomes the 
dirty boson fixed point when Et- = 1. We point out, however, that er = 1 may be quite special. Thus, 
although the qualitative renormalization group flow picture the double epsilon expansion technique 
provides is quite compelling, one should remain wary of applying it quantitatively to the dirty boson 
problem. 

PACS numbers: 64.60.Fr, 67.40.-w, 72.15.Rn, 74.78.-w, 



I. INTRODUCTION 

Quantum phase transitions at zero temperature are 
driven entirely by quantum fluctuations in the ground 
state wavefunction. In many cases a crucial require- 
ment is the presence of quenched disorder. Examples in- 
clude random magnets with various kinds of site, bond^i 
or field^ disorder; the transitions between plateaux 
in the two-dimensional quantized Hall effects?^ metal- 
insulator,'* metal-superconductor, and superconductor- 
insulator^ transitions in disordered electronic systems; 
and the onset of superfluidity of ^He in porous mediajS 

In this article, we will study the quantum transi- 
tions between the superfluid (SF), the localized Bose 



glass (BG), and Mott insulating (MI) phases, motivated 
mainly by the problem of superfluidity of '^He in porous 
media. However, as is often the case in the study of criti- 
cal phenomena, the universality class of this transition, or 
a straightforward generalization of it, also includes other 
physical phenomena, such as many aspects of quantum 
magnetism and superconductivity. 

There have been a number of approaches to study- 
ing the SF~-BG transition: mean fleld theories,^ strong 
coupling)^ large N ^ real space methodSfiS quantum 
Monte Carlo calculationS j^^d^'^'^d^ double dimensional- 
ity epsilon-expansionsfi^ii^'i^, real space renormalization 
group in the limit of strong disorder in one dimensionji^ 
renormalization group in fixed dimension,*^ and in 1 -)- e 
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dimensions. However, none of the above methods pro- 
vide a controlled analytical approach to the critical point 
for d> 1. 

A satisfactory dimensionality expansion about the up- 
per critical dimension for the superfluid to Bose glass 
transition (analogous to the epsilon expansion about 
d = 4 for classical spin problems) does not appear to 
exist. In particular, we previously found that the ap- 
proach based on a simultaneous expansion in both the 
dimension, e,-, of imaginary time (physically equal to 
unity), and the deviation, e ~ A — D, of the total space- 
time dimensionality, D = d + Cr, from four^^ does not 
yield a perturbatively accessible renormalization group 
fixed point, and therefore does not produce a systematic 
expansion for the dirty boson problem. Nevertheless it 
does give one a fairly detailed picture of the renormal- 
ization group fixed point and the flow structures, and 
produce an uncontrolled expansion for the dirty boson 
problem, though further work needs to be done in order 
to understand the analytic structure of the theory as a 
hmction of e,-. Therefore, in spite of its poor convergence, 
the double dimensionality expansion still appears to pro- 
vide the most flexible analytic approach to the study of 
the SF-BG transition, including insight into the symme- 
tries of the fixed point. As detailed below, attention to 
the role of particle-hole symmetry, especially, is essential 
to generating correct RG flow structures. In three (and 
higher) dimensions, the double dimensionality expansion 
is the only method that has been able, with partial suc- 
cess, to access the critical fixed point and obtain critical 
exponentsJi 



A. Particle-hole symmetry 

As stated, it will transpire that an essential ingredient 
that is necessary in order to correctly understand the 
physics of the SF-BG transition is an extra "hidden" 
symmetry, which we call particle-hole symmetry, that is 
present at the critical point, but not necessarily away 
from it. To make this notion precise, we compare the 
following two lattice models of superfluidity: the first is 
the usual lattice boson Hamiltonian, 

where Jij = Jji is the hopping matrix element between 
sites i and j, which we will allow to have a random com- 
ponent; /i is the chemical potential whose zero we fix by 
choosing the diagonal components, Ju, of in such a 
way that Jij = for each i; Si is a random site en- 
ergy with mean zero; Vij = Vji is the pair interaction 
potential, assumed for simplicity to be nonrandom and 
translation invariant; the only nonzero commutation re- 



lations are [ai,aj] = 5ij, and fii = ajoi is the number 
operator at site i. 

The second model is the Josephson junction array 
Hamiltonian, 

Hj = - Jtj cos{4)i - 0j) -I- y^(£t - il)hi 

+ \ YUijninj, (1.2) 

with analogous parameters, but now the commutation 
relations = i5ij. These two Hamiltonians are, in 

fact, very closely related. It is easy to check that if Aq is 
any positive integer then 

a\ = (7Vo + n.)^e"^- 
a, = e-*^'(A^o + "0^ (1.3) 

satisfy the correct Bose commutation relations, and we 
identify hi = Nq + hi. Note, however, that the com- 
mutation relations between (j)i and hi permit hi to have 
any integer eigenvalue, positive or negative, whereas the 
eigenvalues of hi must be non-negative. Therefore it is 
only when Nq is large, and the fluctuations in hi are 
small compared to A'o, that Tij and Tis may be com- 
pared quantitatively: inside the hopping term we may 

approximate a\ k, Nq e'"^% Ui « NqB^'^'^' and make the 
identifications 

Jij = NoJij-; Uij = Vij; e^^£i\ /i = /i - Nq% -I- ivb, 

(1.4) 

where Vq — Va and Vq = i ^^'^ there exists an 

overall additive constant term EqN = {^NqVq — ^Vq — 
Ii)NqN, where A^ is the number of lattice sites. 

Despite this asymptotic equivalence at large Nq, the 
Josephson Hamiltonian (|1.2p has an exact discrete sym- 
metry which the boson Hamiltonian lacks. Thus the con- 
stant shift, h[ = hi + no, where no is any integer, has 
no effect on the commutation relations or the eigenvalue 
spectrum of the hi. The Hamiltonian correspondingly 
transforms as 

nj{h'i} = Hj{h,} + no&o X! + Ne°{nQ,fl), (1.5) 

i 

where e°(no,/i) = no(|noC/o - p.), and Uq = J^jUij. 
The free energy density, fj = In [tre"'^'^-'] , where 
/3 = I/UbT, transforms as 

/j(M) = /j(M-«of/o)+e°("o,M) (1.6) 

independent of the Jij and ii. This implies that the 
only effect of a shift noC/o in the chemical potential is a 
trivial additive term in the free energy which is linear in 
/2. This term serves only to increase the overall density, 
n = by no but otherwise has no effect whatsoever 
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on the phase diagram, which therefore will be precisely 
periodic in //, with period Uq. 

Consider next the transformation = —fii, (p'^ = —(pi- 
The Hamiltonian transforms as 

'>ij[ni,(l)[,£i -p\= Hj[ni,$i, -{ii - /i)], (1.7) 

so that 

fJ{^^Ae^})^fJ{-^^■.{-e^}). (1.8) 

Combining the two symmetries (|1.6[) and (|1.8p we see 
that if all = 0, then for jl = fxi, = ■^kUo, where k is any 
integer, the Hamiltonian possesses a special particle-hole 
symmetry, namely invariance under the transformation 
n'^ = k — fii, (p' = —(/)■ At /i = /ife the density is precisely 
^k per site, and the thermodynamics is symmetric under 
addition and removal of particles (the removal of particles 
being synonymous with the addition of holes). If the Si 
are nonzero, but have a symmetric probability distribu- 
tion, p{£i} = p{—£i}, then the exact particle-hole sym- 
metry is lost, but there still exists a statistical particle- 
hole symmetry at the same special values jlk of fl: self 
averaging will ensure that /j(/ifc, {Ej}) = /j(/ifc, {-Ei})- 
The lattice boson Hamiltonian (jl.ip clearly can never 
possess either form of particle-hole symmetry since the 
hopping term mixes the number and phase in an inextri- 
cable fashion. 



B. Phase diagrams 

In Fig. [T] we sketch zero temperature phase diagrams, 
with and without various types of disorder, in the fl-Jo 
plane, where Jo is a measure of the overall strength of the 
hopping matrix (e.g., Jq = ^i^ij "^y)' simplest 
case of onsite repulsion onlyi^i Uij = UoSij. This phase 
diagram has been discussed in detail previously^^ for 
the lattice boson Hamiltonian, Tis- Here we emphasize 
the features unique to Ti.j, namely the periodicity in jl, 
and the special points (Ik corresponding to local extrema 
in the phase boundaries. 



1. Phase diagram for the pure system 

In Fig. [2a) we show the phase diagram in the ab- 
sence of all disorder. For J^- = the site occupancies 
are good quantum numbers and each site has precisely 
riQ particles for no — ^ < fi/Uo < riQ + ^. The points 
ji/Uo = k + ^ for integer k are 2^ fold degenerate with 
either fc or fc + 1 particles placed independently on each 
site. For > communication between sites occurs 
and the effective wavefunction for each particle spreads 
to neighboring sites (see Fig. (2). We denote by £,{Jo) (to 
be defined carefully later) the range of this spread. One 
can show within perturbation theory,^ however, that for 



sufficiently small, short ranged Jy , there is a finite en- 
ergy gap for addition of particles, and the overall density 
remains fixed at hq for a finite range of fl. 

Consider first fl ^ 0. Then only at a critical value 
Jo.cifl) of Jo does the system favor adding extra particles 
(or holes). Equivalently, for a given Jo there is an interval 
/i-(Jo) < fl < /i+(Jo) of fixed density ng. These extra 
particles may be thought of as a dilute Bose fluid moving 
atop the essentially inert background density, no (see Fig. 
[2]). The physics is identical to that of a dilute Bose gas 
in the continuum, and is well described by the Bogoli- 
ubov modeli^ From this one concludes that the system 
immediately becomes supcrfluid (recall that we assume 
T = 0) with a supcrfluid density ps ^ n — no ^ Jq — Jo,c, 
and an order parameter ipo = {e^^') ^ (n — no)' ^ 
(J — Jo.c)'- The characteristic length in this phase is 

Co = Jq /[fl-fj,±{Jo)]^ ~ (n-no) 2 ^ (Jq_ Jq^^) 2 and 
represents the distance between "uncondensed" particles: 
n — no — l^oP Co"''- This zero-temperature supcrfluid 
onset transition is therefore trivial, in the sense that all 
exponents are mean-field-like. In fact, historically this 
onset was never really viewed as an example of a phase 
transition. 

Furthermore, although all quantities vary continuously 
as Jo decreases toward Jo,c, the actual onset point is en- 
tirely noncritical. Thus, for a given value of Jo within 
a Mott lobe, the interval /i_(Jo) < /2 < /i+(Jo) rep- 
resents a single unique (incompressible) thermodynamic 
state. Incompressibility implies that for the given (inte- 
ger) density, the value of the chemical potential is am- 
biguous. One might just as well set fl — uqUq, its value 
at the center of the lobe. The correlation length, ^(/i, Jq), 
is independent of fl, and remains perfectly finite in the 
Mott phase at Jo,c(/i)- Iii this sense the transition has 
some elements of a first order phase transition. 

The more important transition is the one occurring 
at fixed density, n = no, at fl — uqUq through the tip 
of the Mott lobe at Jo,c(0). At this transition ^(Jq) ^ 
( Jq.c — Jo)"'''""" diverges continuously with a characteris- 
tic exponent, impure- One may show (see Ref.|6iand below) 
that the transition is precisely in the universality class of 
the classical {d + l)-dimensional XY-model. What dis- 
tinguishes this transition from the previous ones is pre- 
cisely particle-hole symmetry: superfluidity is achieved 
not by adding a small density of particles or holes atop 
the inert background, but by the buildup of supcrfluid 
fluctuations within the background, to the point where 
particles and holes simultaneously overcome the potential 
barrier Uq and hop coherently without resistance. The 
exponent impure exhibits itself in the phase diagram as 
well: a scaling argument shows that the shape of the Mott 
lobe is singular near its tip,^ M±('^o) ^ ^{Jo.c — Jo)""""- 
In c? = 2, ^'purc — §, while in d > 3, impure — \- This will 
be discussed in a more general scaling context in Sec. 

HYll 

We have already observed that the lattice boson Hamil- 
tonian (|l.ip never has an exact particle-hole symmetry. 
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FIG. 1: (Color online) Schematic zero temperature phase diagram for the Josephson junction model, (a) The model without 
disorder, showing the periodic sequence of Mott lobes (MI) and the superfluid phase (SF). The transitions to superfluidity at 
the tips of the Mott lobes, J = , are special and are in the universality class of the {d + l)-dimensional XY-model. The 
points fJ.k/Uo ~ k/2 have an exact particle-hole symmetry. As described in the text, there is a corresponding singularity, 
determined by the correlation length exponent i^^Ji, in the shape of the lobe near its tip. (b) The model with site energy 
disorder whose distribution is supported on the interval — < Ei < A^, showing the now shrunken (or even absent, if the 
disorder is sufficiently strong, 5 = Ae/Uo > |) Mott lobes, and the new Bose glass phase that now intervenes between them 
and the superfluid phase. As indicated, the boundaries of the Mott lobes are determined entirely by 5 and their pure system 
counterparts. The transition to superfluidity now takes place only from the Bose glass phase (BG), and it is believed that 
the nature of this transition is independent of where it occurs, even at the special points jj,k which now have only a statistical 
particle-hole symmetry, (c) Phase diagram for the model with only nearest neighbor bond disorder, Jij = Jo(l -I- AJtj) whose 
distribution has support —5- < A Jij < 5+, with 5_ < 1, showing again shrunken Mott lobes (which will vanish entirely if 
A Jij is unbounded from above). The boundaries of the Mott lobes are now determined entirely by S+ and their pure system 
counterparts. At integer filling the Bose Glass phase turns into an incompressible random rod glass (RRG). The transition is 
in the universality class of the (d -I- 1) -dimensional classical random rod problem. Away from integer filling the compressibility 
turns on continuously, and the Bose glass phase reappears. The transition is once more in the universality class of the generic 
Bose glass to superfluid transition. The random rod temporal correlation length exponent Vrfi ~ zqUq (expected to be greater 
than unit in d = 3^^) now describes the shape of critical line near the commensurate point. In a neighborhood of half fllling 
(half-integer fi/Uo), whose size must depend on the precise distribution of A Jij, for d > 1,^^ the system remains superfluid all 
the way down to Jo = 0: in the absence of site disorder the the exact particle-hole symmetry that is restored at half fllling 
guarantees delocalized excitations. 



Nevertheless, one still has Mott lobes (now asymmetric 
and decreasing in size with increasing no) for each integer 
density, and a unique extremal point, [Jo,cino), fidno)], 
at which one exits the Mott lobe at fixed density n = n^. 
One may show'' that the transition through these ex- 
tremal points is still in the [d + l)-dimcnsional XY uni- 
versality class, and that particle-hole symmetry must 
therefore be asymptotically restored at the critical point. 
The difference now is that there is a nontrivial balance 
between the densities of particle and hole excitations, and 
the interactions between them. The position of the criti- 
cal point is no longer fixed by an explicit symmetry, but 
must be located by carefully tuning both the hopping 



parameter and the chemical potential. 

This phenomenon of "asymptotic symmetry restora- 
tion" at a critical point is actually fairly common (and 
we shall encounter it again below). For example, though 
the usual Ising model of magnetism has an up-down spin 
symmetry, the usual liquid-vapor or binary liquid crit- 
ical points do not. However the Ising model correctly 
describes the universality class of the transition, and one 
concludes that the up-down symmetry must be restored 
near the critical point. Similarly, the p-state clock model 
with Hamiltonian 



Ti. = — cos 



27r ■ 

—{q^-<lJ) 



= l,2,...,p, (1.9) 
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FIG. 2: (Color online) Schematic illustration of lattice bosons 
near unit filling n = 1 at some value of the hopping amplitude 
in the interval < Jo < Jo,c- In the Mott insulating phase 
(n = 1), the effective wavefunction of each particle spreads 
only a finite distance ^(Jo), and the state is insulating. For 
n > 1 (n < 1), the extra particles (holes) travel freely within 
the essentially inert background, and the state is superfluid 
for arbitrarily small |n — I|. 



which may be thought of as a kind of discrete XY-mode\, 
has for sufficiently large p (specifically, p > 4 in d — 2; 
clearly p = 2 corresponds to the Ising model and p — 3 
to the three-state Potts model) a transition precisely in 
the XF-model universality class. Note, however, that in 
the ordered phase, corresponding to the zero temperature 
fixed point, the order parameter will (for d > 2) sponta- 
neously align along one of the p equivalent directions, qi, 
breaking the Xy-symmetry and generating a mass for 
the spin-wave spectrum (more interestingly, in d = 2 a 
power-law ordered Kosterlitz-Thouless phase exists for a 
finite temperature interval below the transition, with a 
second transition to a long-range ordered phase taking 
place only at a lower temperature^!) . This latter prop- 
erty is not relevant in the present case since breaking 
particle-hole symmetry does not break the symmetry of 
the order parameter: the nature of the superfluid phase 
is unaffected. 



2. Phase diagrams with disorder 

We will consider two types of disorder: (i) onsite disor- 
der in the ii, and (ii) disorder in the hopping parameters, 
Jij. li pL ^ jlk for any k, that is if particle-hole symme- 
try is broken, we expect the two types of disorder to 
yield the same type of phase transition. In renormaliza- 
tion group language, each in isolation will generate the 
other under renormalization. If fl = flk and the ii have 
a symmetric distribution about zero so that the Hamil- 
tonian possesses a statistical particle-hole symmetry, the 
obvious question is whether or not the transition in this 
case is different from the one in the presence of generic 
nonsymmetric disorder. We shall argue below that it is 
not, i.e., that breaking particle-hole symmetry locally is 
not substantially different from breaking it globally, and 
that in fact statistical particle-hole symmetry is asymp- 




FIG. 3: (Color online) Schematic illustration of lattice bosons 
near unit filling in the presence of bounded disorder. If the 
disorder is not too strong (A^ < Uo/2), there is still a Mott 
phase with a finite energy gap for adding or removing par- 
ticles. Unlike in the pure case (Fig. [2]), superfiuidity is not 
generated immediately with the addition of particles. For 
sufficiently small n — 1, the additional particles are Ander- 
son localized by the residual random background potential 
of the effectively inert layer. The finite compressibility dis- 
tinguishes this Bose glass phase from the Mott phase. The 
superfluid critical point /iaf(Jo) occurs only once the added 
particles have sufficiently smoothed the background potential 
that its lowest lying states become extended. 



totically restored at the critical point Only \i jl = fl^ 
and £i = does the disorder fully respect particle-hole 
symmetry. We shall see that in this case the transition is 
entirely different, lying in the same universality class as 
the classical [d + l)-dimensional XF-model with colum- 
nar bond disorder, precisely the kind of system addressed 
in Ref . 3 

In Fig. [ijb) we sketch the phase diagram in the pres- 
ence of site disorder, whose distribution is supported on 
the finite interval — < < A^. We see that the Mott 
lobes have shrunk (and may in fact disappear altogether 
for sufficiently strong disorder), and a new Bose glass 
phase separates these lobes from the superfluid phase. ^ In 
this new phase the compressibility is finite, but the parti- 
cles do not hop large distances due to localization effects: 
particles on top of the inert background still see a resid- 
ual random potential, whose lowest energy states will be 
localized (Fig. [31). As particles are added to the system, 
bosons will tend to fill these states until the residual ran- 
dom potential has been smoothed out sufficiently that ex- 
tended states can form, finally producing superfluidity?^ 
As argued above, the nature of the superfluid transition 
is the same everywhere along transition line. 

As indicated in the figure, the boundaries of the 
Mott lobes are determined entirely by the pure system, 
together with A^.^ The upper half of the boundary, 
/24.(Jo,Ae) = /i+(Jo,0) - Ae, is pushed down by A^, 
while the lower half, /i+(Jo, Ag) = /i_(Jo,0) -I- A^, is 
pushed up by A^. This result, which relies on the exis- 



tence of large, rare regions of nearly uniform superfluid, 
will be derived in Sec. IIIII 

Finally, in Fig. [Ijc) we sketch the phase diagram in 
the presence of bond disorder. For simplicity we consider 
here only nearest neighbor hopping parameterized in the 
form Jij = Jo{l + AJij), with [AJy]^^ — 0, whose distri- 
bution is supported on a finite interval —S- < A Jij < 6+, 
with (5_ < 1. The Mott lobes have again shrunk (and 
may in fact disappear altogether for unbounded disor- 
der, 6+ oo), and glassy phases again separate these 
lobes from the superfluid phase. At non-integer density, 
the glassy phase is the compressible Bose glass. How- 
ever, the existence of an exact particle-hole symmetry at 
integer density changes the nature of the glassy phase 
there, turning it into an incompressible random rod glass 
(RRG). In both cases, localization effects destroy super- 
fluidity over a finite interval.— As argued above, at non- 
integer density, the universality class of the BG-SF tran- 
sition is of the same as for the site disorder model, while 
at integer density the RRG-SF transition is in the special 
random rod universality class. 

The random rod temporal correlation length exponent 
I'r.o = zoVq exhibits itself in the phase diagram: as also 
discussed in Sec. IIVBI the superfluid transition line has 
a singularity, ^±(Jo) ~ M^f^^ - JoT""" , in the neigh- 
borhood of the particle-hole symmetric points (note that 
Zpure = !)■ One expects i^r.o > 1 in d = 3, yielding the 
pictured cusps. 

The boundaries of the Mott lobes are again deter- 
mined entirely by the pure system, together with (5+ . The 
boundary, Jo,c(Aj'^+) = •-'o,c(Mi 0)/(l + S+), is this time 
scaled to the left by Aj. The result again relies on the 
existence of large, rare regions of nearly uniform super- 
fluid, and will also be derived in Sec. IIIII 

At half fining, where /x = /ifc = (A; -f ^)Uo is a half- 
integer, an exact particle-hole symmetry is restored, and, 
for d > 1,— the superfluid phase can penetrate all the 
way to zero hopping. In general the superfluid phase will 
survive on a finite interval in density (which maps into 
the single point /2 = /ifc at Jo = 0) around half filling, 
whose size depends on the precise distribution of Jij.^^ 



C. Criticality and restoration of statistical 
particle-hole symmetry 

In Sec. [TT] various functional integral representations 
of the Hamiltonians (jl.ip and (|1.2p will be introduced. 
In order to discuss, in the most transparent fashion, the 
role of various symmetries at the superfluid transition, 
we summarize in Table U the basic classical continuum 
-0* models that may be abstracted from these represen- 
tations. The phase of ?/; = \ip\e'^'^ represents the Joseph- 
son phases in (jl.2p . The control parameter represents 
the hopping strength Jq, while represents the chemi- 
cal potential /2. Quenched hopping and site energy dis- 
order are correspondingly represented, respectively, by 
the random fields (5r(x) and (5g(x). The fact that they 
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FIG. 4: (Color online) Schematic illustration of random rod- 
like structure generated by quenched disorder in quantum 
models. The cylinders represent, for example, a picture of 
random regions of enhanced or suppressed hopping strength. 



are r- independent generates rod- like disorder (Fig. d]). 
These models are intended to be used in the vicinity of 
the n = Mott lobe, hence in the neighborhood of 
zero, since fluctuations in the amplitude destroy any 
translation symmetry in analogous to (|1.6p . and the 
nonzero lobes are no longer symmetric. 

The Lagrangian Cq = Ci{gf) = 0) generates the clas- 
sical {d + l)-dimensional XY critical behavior along the 
line /i in Fig. [IJa): the superfluid transition occurs 
for decreasing tq at a critical value ro,c- Nonzero go in >Ci 
generates the remainder of the Mott lobe, r-o,c(5o), corre- 
sponding to the remainder of the phase diagram in Fig. 
[ija). The transition at nonzero go is in the universality 
class of the dilute Bose gas superfluid onset transition 
(described by the Bogoliubov model) in d dimensionsj^i 
The usual coherent state representation of the (contin- 
uum version of the) boson Hamiltonian (jl.ip is obtained 
by dropping the |9r'0P term in Ci and setting g^ = 1. 
This demonstrates explicitly the lack of a simple inter- 
polation between the particle-hole symmetric and asym- 
metric models in the original boson Hamiltonian. 

The Lagrangian £2 = ^^{ga = 0) incorporates hop- 
ping disorder in the form of r-independent disorder 5r in 
ro, and represents the previously analyzed (particle- hole 
symmetric) random rod model.— The model describes 
the transitions along the line g^ = Q m Fig. [ijc) , includ- 
ing the incompressible random rod glass (RRG) separat- 
ing the tip of the Mott lobe from the superfluid phase. 
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Pure PH-sym [{d + l)-dimensional XY model]: 

£o = - / d'^x J dr + lia.V'P + IrolV^p + \uo\^\'} 

Pure PH-asym [(i-dimensional dilute Bose gas]: 
Ci = -J d'^xjdr {i|V^|2 - iridr - 9o?i^ + irol^p + |«o|^r} 

PH-sym RR [{d + l)-dimensional classical random rod model]: 
£2 = - / rf'^x/ {i|V^|2 + l\drij\' + |[ro + 5r(x)]|^|2 + l^ol^r} 

PH-asym RR [{d + 1) -dimensional incommensurate random rod model]: 
Cs = -Jd''xJ rfr {i| - iridr - gof^ + \[ro + '^r(x)]|^|2 + iwolV'l'} 

Statistical PH-sym [commensurate dirty boson problem]: 
£4 = - / d'^x^dr {i|V^|2 - \r\dr - ^<7(x)]2V' + \\r^ + (5r(x)]|V'|2 + \uM\'} 

Generic PH-asym [incommensurate dirty boson problem]: 
£5 = - / c^'^x/ {i|V^/^|2 - \r\dr -g^- 5g{^)f^ + \[ro + 5r(x)]|V;p + iwolV'!'} 

TABLE I: ip'^ representation of models with various types of disorder and various degrees of particle-hole symmetry. The 
coefficients of IVi/ip and |9r'i/'P have been normalized to i. The control parameters ro and go are analogous to Jo and fi, 
respectively. Disorder in the hopping strengths is represented by 5r, while that in the site energies is represented by 5g. Both 
are independent of r. In field theoretic treatments, both are taken as quenched Gaussian random fields with zero mean and 
delta-function correlations characterized by variances and Ag , respectively. Disorder in the other parameters (including the 
unit gradient-squared coefficients) may also be introduced, but produces no new critical behavior. 



Nonzero go in -Ca generates the remainder of the phase 
diagram in Fig. [ijc). We will show in Sec. IIIII that the 
breaking of particle-hole symmetry implies that the RRG 
becomes the compressible Bose glass exists at all nonzero 
go — in the renormalization group sense, nonzero go, to- 
gether with nonzero (5r, generates a finite 5g even if it 
vanishes to begin with. 

Under the condition that the distribution of 6g is sym- 
metric, he Lagrangian £4 = £5 (50 = 0) maintains a sta- 
tistical particle-hole symmetry. It is believed, however, 
that the superfluid transition at zero and nonzero 50 are 
in the same universality class. In renormalization group 
language, go must flow to zero at large scales, so that £4, 
not £5, describes the critical fixed point. This is the tech- 
nical definition of the restoration of statistical particle- 
hole symmetry. We will confirm this picture within the 
double-epsilon expansion in Sec. IVII 

The significance of this symmetry restoration at the 
critical point is the following. In Ref. [l^ only the boson 
Hamiltonian p.ip , and corresponding coherent state La- 
grangian, was considered, and a fixed point sought only 
in the space of parameters accessible to this Hamiltonian. 
Such a fixed point can therefore never possess a statistical 
particle-hole symmetry. However if the true fixed point 
does possess this symmetry, it is clear that it must then 
lie outside the space of boson Hamiltonians of the form 
(II. ip (or Lagrangians without a ISrV'P term). Accessing 
this fixed point requires an enlargement of the parameter 
space so that both particle-hole symmetric and asymmet- 



ric problems may be treated within the same model. The 
Josephson junction array Hamiltonian, (|1.2p . from which 
all of the Lagrangians listed in Table |T] may be derived, 
satisfies this requirement. 

Using £4 and £5, we shall see in Sec. |VT] that, in- 
deed, a new statistically particle-hole symmetric dirty 
boson fixed point can be identified, and that all technical 
problems encountered in Ref. [l^ may then be avoided. 
Specifically, certain diagrams that were ignored in Ref. 
Ilg are in fact important and accomplish the required en- 
largement of the parameter space. There is, however, a 
price: this new fixed point is not perturbatively accessible 
within the double epsilon expansion. Thus, the random 
rod problem may be analyzed for small eind = 4— e — e,- 
dimensions if the dimension, of time is also smallj^ 
We find, however, that for small e and e,-, site disor- 
der is irrelevant at the random rod fixed point, and that 
full particle-hole symmetry is therefore restored on large 
scales close to criticality. This remains true for suffi- 
ciently small tr < e-%{D), with e%{D) = 0(1). Only for 
Et- > e%{D) does a new fixed point appear which breaks 
full particle-hole symmetry. Extrapolation of the critical 
behavior associated with this new fixed point, though un- 
controlled, shows significant similarities to some features 
of the known behavior at e^. = 1. To lowest nontrivial 
order in we find t%{D = 4) ~ ^ (D = 4, hence e = 0, 
corresponding to d = 3 at e^. = 1). This value of is 
significantly less than unity, and leads us to hope that es- 
timates based on these extrapolations from small Et- are 
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not too unreasonable. 



D. Outline 

The outline of the remainder of this paper is as follows. 
In Sec.|TT]we introduce various useful functional integral 
formulations for the thermodynamics of the Hamiltoni- 
ans (II. ip and (|1.2p . We begin in Sec. Ulllby considering 
the role of particle-hole symmetry in the nature of the 
excitation spectra of the glassy phases. Using a phe- 
nomenological model in which we view the structure of 
the random rod and Bose glass phases as a set of ran- 
dom sized, randomly placed, isolated superfluid droplets, 
we focus on the density and compressibility and exam- 
ine how they vanish as full particle-hole symmetry is re- 
stored. The droplet model also confirms the relation be- 
tween the boundaries of the pure and disordered Mott 
lobes in Figs. [TJ In Sec. |TV] we begin focusing on the 
critical point through various phenomenological scaling 
arguments. In particular we identify a new crossover ex- 
ponent that describes the relevance of particle-hole sym- 
metry breaking perturbations to the random rod critical 
behavior. We revisit the original arguments-iSS. for the 
dynamical exponent scaling equality z — d, showing that 
they can be violated, and hence that z may remain an 
independent exponent This is consistent with recent 
quantum Monte Carlo simulations in d = 2 that find 
z — 1.40±0.02.1'* We also discuss the asymptotic restora- 
tion of statistical particle-hole symmetry at the dirty bo- 
son critical point. In Sec. |V| we illustrate all of these 
ideas using an exactly soluble one-dimensional model. 
The analysis is very similar to that of the Kosterlitz- 
Thouless transition in the classical two-dimensional XY- 
model. In Sec. IVII we generalize the previous analyses to 
general Et- 7^ 1, observing along the way some apparent 
pathologies that make e,- = 1 very special, leading one 
to question how smooth the limit e,- — > 1 might be. For 
example, the Bose glass phase has finite compressibility 
for Ct = 1, but is incompressible for all e,- < 1- It is dis- 
tinguished from the Mott phase only by having a diver- 
gent order parameter susceptibility, Xs- We then intro- 
duce the Dorogovtsev-Cardy-Boyanovsky double epsilon 
expansion formalism^^ and derive the results outlined in 
the previous subsection. Finally, two appendices outline 
the derivations of various path integral formulations and 
duality transformations used in the body of the paper. 



II. FUNCTIONAL INTEGRAL 
FORMULATIONS 

In order to obtain a formulation of the problem more 
amenable to analytic treatment, we turn to functional 
integral representations of the partition function. It will 
turn out to be important to have an exact representa- 
tion. Representations which involve dividing the Hamil- 
tonian into two pieces, Ti. = Tio + Tii, then using the 



Kac-Hubbard-Stratanovich transformation to decouple 
Til, generate effective classical actions with an infinite 
number of terms, which must be truncated at some finite 
order.® In addition, such representations work only when 
fl lies within a Mott phase when Jg — 0, and hence break 
down for unbounded, e.g., Gaussian, distributions of site 
energies. We turn instead to representations obtained 
from the Trotter decomposition (see Add. E]). 



A. Lattice boson model 

For the lattice boson model, the coherent state repre- 
sentation is most appropriate, and yields a classical La- 
grangian 



Cb = 



dr 



^*{T)drMr)~nB{i^:{T),Mr)} 



(2.1) 

where 'H{ip*{T),ipi{T)} is obtained by substituting the 
classical complex variable ^'ti'^) for the boson site anni- 
hilation operator [and ^1{t) for the creation operator 
a|] wherever it appears in the (normally ordered form 
of the) quantum Hamiltonian. The partition function 
IS given hy Z = tr'^[e^], where tr'^[-] is an unrestricted 
integral over all complex fields, ipiir). Notice that the 
only term that couples different time slices is the "Berry 
phase" il)*dr4' term which arises from the overlap of two 
coherent states at neighboring times. This should be 
contrasted with the spatial coupling, i j (es- 
sentially a discrete version of ?/'*V^^), which appears in 
TLb- The imaginary time dimension is therefore highly 
anisotropic. This anisotropy is increased further if disor- 
der is present since the £i and Jij are r-independent: the 
disorder appears in perfectly correlated columns, rather 
than as point-like defects, in [d + l)-dimensional space- 
time. 

If the '4>*dri' term were replaced by ip*d'^ip, only the 
disorder would contribute to the anisotropy (the fact that 
the coefficients of ip*d'^ip and iJ;*V'^ip are different is not 
important, and may be cured by a simple rescaling). The 
model becomes precisely a special case in the family of 
classical models with rod-like disorder treated in Ref. [isl . 
The linear time derivative term in (|2.ip is more singular 
than a term with a second derivative in time. It should 
therefore not be too surprising that its presence leads to 
new critical behavior. 

Yet another crucial property of the il}*dr4' term is that 
it is purely imaginary: 



^0 



drip* drip 



dT'ip*dT'4'j 



(2.2) 



where integration by parts and periodic boundary con- 
ditions have been used. Therefore the statistical factor, 
e^^ , used to compute the thermodynamics is in general 
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a complex number, and leads to interference between dif- 
ferent configurations of the V'i (''')■ Unlike that with cou- 
pling ip*d'^ip, the resulting model therefore does not cor- 
respond to any classical model with a well defined real 
Hamiltonian in one higher dimension. In fact, it is pre- 
cisely this property that reflects the particle-hole asym- 
metry in the model. Interchanging particles and holes is 
equivalent to interchanging i^lir) and ipi^r). The ip*dr'ip 
term changes sign under this operation, while TCb[iI'* t'4'] 
is unaffected. Thus although £b is invariant under the 
combination (known as time reversal) of complex con- 
jugation and r — r, the boson model always violates 
each separately. 



B. Josephson model 

Consider now the canonical coordinate Lagrangian for 
the Josephson array model (see App. lAlfor a derivation): 



J X] ^'J cos[0,(t) - (t)j{T)] 



i(t>j{T) 



(2.3) 



with partition function Z = tr'^e'''. Notice that the 
linear time derivative, 4>i, now appears in a much more 
symmetric looking fashion. If // — = 0, Cj is particle- 
hole symmetric and real: 



£j[M + e, =0] 



dr 



^ JyC0s[(/),(T) - (/)j(r)] 

id 



(2.4) 



and takes precisely the form of a classical XF-model in 
(d + 1) dimensions. 

The periodicity, ()1.6p . of the phase diagram is a con- 
sequence of the periodicity of the 0^ in r: substituting 
/i — uoUq for p, and multiplying out the {U~^)ij term 
yields 

-Cji/i - no?7o] = Cj[p]-inQ / dr'S^'PiiT) 

Jo 

+ /3iVeO(no,M). (2.5) 

However dr^iir) = 2-Kmi and e'^''™*"" = 1, so the 
second term simply drops out of the statistical factor, 
e'-'-' , and we recover the free energy identity (|1.6p . Notice 
that if jl — f-i = we obtain a statistical factor 



aCj[fi-ei = ^Oo\ 



(-1) 



£j[/i-ei=0]+/3Are''(i,0) 



, (2.6) 



which, though real, is not always positive. Although this 
Lagrangian is also particle-hole symmetric, it too does 
not correspond to the Hamiltonian of any classical model. 
This model is very different from that with jl~ei = 0. For 
example, as shown in Fig. [H^c), it always has superfluid 
order at T = 0, for arbitrarily small Jy, as opposed to 



2.4|) which orders only for sufficiently large Ji^ 



C. Josephson model for general 

For later reference, note that, in contrast to (|2.1[) . the 
Lagrangian (j2.3p has an obvious generalization to nonin- 
teger dimensions of timei^ If Cr is the dimension of time, 
we simply write 



= j^^ d'-r jlj; 4- cosiM-r) - 0,(t)] + ^ Jli^-'h [*V.</),(t) + - e,] • [*V.(/.,(t) + M - Cj] | , (2.7) 



where t, fj. and are e^-dimensional vectors: r = (ri, . . . , r^^), etc. 

Renormalization group calculations are performed most conveniently on Lagrangians, such as (|2.ip . which are 
polynomials in unbounded, continuous fields and their gradients. Therefore wc would like to convert (j2.7p to such a 
model, while retaining the essential physics. If we write '4'i{T) — 6"^'*^"^^ then (|2.7p may be written 



'-J 



[ ^'^^ ^J^3W^{T)iJ,{T) + C.C.] + \ Y^{U-%i^*{T){V r + /X - e,)^.(r) • V*(t)(V. + ^l - e,)ij,{T)^ . 



For onsite interactions only, J7y = J7o(5.y, the second term which is conveniently quadratic in V'- We now relax 
simplifies to the assumption jV'il = 1, employing instead the usual 
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r|-0p + Landau-Ginzburg- Wilson weighting factor, 
obtaining finally 



Y^[rmT)?+vmT)\' 



(2.10) 



This model retains the exact particle-hole symmetry at 
/X + Ci = 0, but loses the precise periodicity of the phase 
diagram when Et- = 1: thus the second term in (|2.5p now 
becomes 



no 



rl3 

Jo 



(2.11) 



[compare the first term in (j2.ip ] which reduces to the 
previous form if iT/^il = 1. However if [tpil fluctuates, as 
in (j2.10p , this term is no longer a perfect time derivative 
and will not integrate to a simple integer result. We will 
therefore only use (|2.10p near fl ~ when we study the 
role of particle-hole symmetry near the phase transition. 



D. Continuum models 

The field theoretic LGW-type Lagrangians listed in Ta- 
ble [T] follow (for er — 1, but with obvious generalizations 
to general er) from the continuum limit of (|2.10p (and 
its various special cases), in which the nearest neigh- 
bor hopping term maps to jV^/ip. Space and time are 
also rescaled to produce unit coefficients of the 9tV'P 
and IVV'P terms. As a result, the hopping disorder now 
maps to disorder in the \ip\'^ coefficient. As is standard, 
the mapping is not exact, but is intended only to produce 
a minimal model that preserves the basic symmetries of 
the original, so that its phase transitions lie in the correct 
universality class. The control parameters and have 
only a rough correspondence with Jq and pt, but never- 
theless generate phase diagrams with the same topology. 



III. PARTICLE-HOLE SYMMETRY AND THE 
EXCITATION SPECTRUM OF THE GLASSY 
PHASES 

In this section we will consider the nature of the non- 
superfluid phases in the presence of the two types of dis- 
order, and Jij (to simplify the notation we henceforth 
drop the tildes on the Josephson junction model param- 
eters). Recall that for the boson problem the produce 
a Bose glass phase,— with a finite compressibility and a 
finite density of excitation states at zero energy. We shall 
contrast this with the case of the particle-hole symmetric, 
random Jy model, which we will show has a vanishing 



compressibility and an excitation spectrum with an ex- 
ponentially small density of states, p{e) ~ e""^"/^, i.e. a 
"soft gap." We will show that the compressibility is pre- 
cisely the spin-wave stiffness in the time direction, which 
therefore vanishes in the "symmetric glass," but is finite 
in the Bose glass. This yields an upper bound zq < d 
for the dynamical exponent at the particle-hole symmet- 
ric transition. An effective lower bound on zq may be 
obtained by demanding that particle-hole asymmetry be 
a relevant operator at the symmetric transition. This 
is a necessary condition in order that the particle-hole 
asymmetric transition be in a different universality class 
from the symmetric one. We will obtain estimates for 
this lower bound within the double e-expansion in Sec. 

m 



A. Superfluid densities or helicity modulii 

We begin by defining the superfluid density — the "he- 
licity modulus," or "spin wave stiffness," in classical spin 
models. This quantity is computed from the change in 
free energy under a change in boundary conditions. Con- 
sider a box-shaped system with sides La, a = 1, . . . , Z?. 
We are interested m D ^ d+1 and Ljj = p. We say that 
ip obeys 0Q-boundary conditions if 

il:{xi, . . . ,Xa + La, ■ . ■ ,xd) = e''^ i>{xi, . . . , Xa, ■ ■ ■ , X d) 
i>{xi,. . . ,X,3 + L,3, . . . ,Xd) = 'iJj{xi,...,Xi3,...,Xd), 

P^a (3.1) 

i.e., a twist angle 9 is imposed in the a direction, while 
periodic boundary conditions are maintained in all other 
directions. Let 



Intr 

^ Vd 



,Vd^\{ Lp, 



(3.2) 



13=1 



where is the Lagrangian, be the free energy obtained 
using 0Q-boundary conditions. We may define 



7^(a;i,...,xc) = e-^'""/^°V(^i 



(3.3) 



which obeys periodic boundary conditions in all direc- 
tions. In all cases of interest one may write 



C'-[iP]^C''[i^]+5C[t,eiLa] 



(3.4) 



where 5C may be expanded as a Taylor series in powers 
of 9 / La, and superscript "0" denotes periodic boundary 
conditions. Thus 



1 



{5C) + \{5C'') 



, (3.5) 



where {6£'^)c — {SC^) — {SC)"^ , and the averages are with 
respect to >C°[i/']- Equation ()3.5p yields a series of terms 
in powers of 9/ La, and we use the notation 



= .l^p^ + -(— 

La 2 \La 



(3.6) 
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to define the coefScients in this series. This only makes 
sense for \0\ < tt, since it is clear from the definition that 
is periodic in 0^ with period 27r. We shall see that the 
first term, which in most previous cases was completely 
absent, arises from particle-hole asymmetry. The second 
term defines the helicity modulus, Tq, in the direction a. 



1. Boson model helicity modulus 

Let us now turn to specific cases with Lagrangians de- 
fined by (12. ip and (|2.3|) . In both cases, when a is a spa- 
tial coordinate {a — \, . . . ,d) the sensitivity to boundary 
conditions comes only from the hopping term so that for 
the boson model [see (|l.ip and (|2.ip ]. 



5Ci 



1 

2 7o 



2^0 Jq 

q2 r-P 

4L~ 



dT^J,,(x«-a;«)[^*^,-c.c.] 



J f drY, J,, {xt - x'^nr.i^, + C.C.] 



(3.7) 



2. Josephson model helicity modulus 

The Josephson Lagrangian yields precisely the same 
expressions if one identifies V'il'r) = e"^*'"^'. Thus, (3.8) 
becomes 



T — - 



X (sin[(/),(r) - (/.j(r)] sin[0fe(O) - MO)]) 
Y,MO^\{cos[MO) ' MO)]) 



a = 1, . . . , d 



(3.10) 



and l|3.9p becomes 



rP p 
T„ = -U^a^ J drY, [(si 



sin[(/>i+x^(T) - (t>i{T)] 



X sin[0x„(O)-</)o(O)]) 



2Ja" 



cos[</)x„(O)-0o(O)]) 
a = 1, . . . , d 



(3.11) 



5. Temporal helicity modulus and compressibility 



which, due to the assumed vanishing of the net current 
under periodic boundary conditions, yields pa = 0, a = 



1, . . . ,d, and 



X ([^*(r)V',(r)-c.c.][V'^(0)Vo(0) 
+ ^E['^'o(^")'(^*(0)^o(0) + c- 



c.) 



a 



(3.8) 



where [-Jav denotes an average over the disorder (we as- 
sume self averaging). For nonrandom J^ , with nearest 
neighbor hopping J only, p.8p reduces to 



t2 2 



drY [{[r^{r)^aMr) ~ Mr)dc^4ir)] 



X [ro{Q)dMo) - i^^io)^o.^:m) 



+ J a' 



, a = \, . . . ,d, (3.9) 



where = '^i+stc ~ V'jj is the lattice spacing, and, 

in an obvious notation, i -I- Xq, labels the nearest neigh- 
bor lattice site in direction a. We recognize this as the 
discrete version of the usual definition of in terms of 
the current-current correlation function. 



Consider next the stiffness in the time direction. We 
will show that it is precisely the compressibility, n = 
— -Q^- To see this, note that only terms with time deriva- 
tives are sensitive to 6'r-boundary conditions. In the Bose 
case, Eq. ()2.ip . we obtain 



6£ = '-^ iyri::ir)Ur) 



(3.12) 



which corresponds precisely to an imaginary shift, fj,' = 
p + i^, in the chemical potential. Similarly, for the 

Josephson case, Eq. (|2.3p . we define (f>i{T) — 4>i{T) ~ ^t, 
leading to exactly the same chemical potential shift. 
Thus in both Cb and Cj the time derivatives appear 
with the chemical potential in just the right way to give 
rise to what amounts to the Josephson relation between 
the time derivative of the phase and changes in the chem- 
ical potential. We immediately conclude that the series 
p.6p takes the form 



i0 df 



ey 52/° 



13 J dfi^ 



(3.13) 



where p = is the number density, and, as promised, 

we identify T^. — k. 

Classical intuition tells us that Tq should be nonzero 
only when the model has long range order in the phase 
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of the order parameter, i.e. only in the superfluid phase. 
Ahhough this statement is true for the spatial directions, 
a = 1, . . . ,d, this is not necessarily true for a = t. Our 
intuition regarding ^He in porous media would lead us 
to be very surprised if the system were incompressible, 
K = 0, throughout the nonsuperfluid phase. Thus there 
should be no barrier to the continuous addition of parti- 
cles to the system, even when it is completely localized 
(only Mott phases, in which disorder is unimportant, are 
incompressible because the density is pinned at special 
values commensurate with the lattice^). The Bose glass 
phase is therefore rather special in that the order param- 
eter has a temporal stiffness, = k > 0, even when 
there is no spatial stiffness, T^; oc ps = 0. Our classical 
intuition breaks down because the Lagrangian is typically 
not real and, as discussed earlier, does not have a proper 
classical interpretation. The droplet model picture de- 
veloped below will make clear the origin of this special 
partial order. 

The particle-hole symmetric model described by (|2.4p . 
however, does have a classical interpretation, and de- 
spite the fact that the Jij are random and the model 
anisotropic (the disorder being fixed in time) it would be 
surprising if the disordered phase possessed long range or- 
der in time. Thus we expect k to vanish when ps does, so 
that the glassy disordered phase is incompressible. This 
is permitted because the particle-hole symmetry now dic- 
tates that the density be an integer. What distinguishes 
this glassy phase from the Mott phase, however, is that 
K is not zero for an entire interval of /i, but vanishes only 
for the special value /i = where particle-hole symmetry 
holds. 



4- Continuum model helicity moduli 



bonds only, with all A Jjj independent random variables 
with zero mean. Let Jc be its critical point, and let 
be the critical point when all 6 Jij — (note that it is 
entirely possible that Jc < J^, since random fluctuations 
can sometimes compete with Mott phase commensura- 
tion effects and thereby enhance superfluid order^^). In 
the latter, nonrandom case, the transition is from a Mott 
insulating phase for Jq < to a superfluid phase for 
Jq > J°. Suppose now that < AJ^ < S+ is bounded 
from above (as well as, trivially, from below) with the 
essential supremum (i.e., the largest value of A Jij achiev- 
able with finite probability density). Then for Jq such 
that Jo + S+ < J°, all Jij are smaller than J°, and the 
system must have a Mott gap — reducing any set of the Jij 
from an initially uniform value in the Mott phase can only 
enhance its stability. However, for Jc > Jq > J°/(l-|-(5+) 
one will form, via probabilistic fluctuations, exponen- 
tially rare, but arbitrarily large regions of bonds in which 
all Jij > Jc- These regions therefore represent finite 
droplets of superfluid — increasing any set of the Jij from 
an initially uniform value within the superfluid phase can 
only enhance the stability of the superfluid phase phase. 
It is here that the r-independence of the Jij becomes 
critical — in the classical interpretation these droplets are 
one-dimensional cylinders (see Fig. [¥]) with arbitrarily 
large cross-section, made of material that would be fer- 
romagnetically ordered in the bulk. The fact that these 
regions are already infinite along one dimension clearly 
enhances magnetic ordering more than would finite (zero- 
dimensional) pieces of ferromagnet. We shall see now 
that these droplets generate Griffiths singularities^ that 
close the Mott gap. It then follows that the Mott phase 
boundary must lie precisely at J°/(l -|- (5+), as shown in 
Fig. He). 



For the continuum models listed in Table [H the simple 
squared gradient term produces in all cases the isotropic 
result 



T — 



(3.14) 



d'^X / dT[(^*5„V(x,T)^*a„^(O,0))]^ 




for the spatial helicity moduli. The temporal phase twist 
response is given by (j3.14p with go replacing fj,: p = 

Clearly, for Cq, £2, and £4, 
after taking the derivatives. 



.9fl 
dgo 

one should set go 



B. Droplet model of the glassy phases 

Let us now understand in detail how these two differ- 
ent behaviors merge with each other in the full phase 
diagram, Fig. [1] and in particular confirm the positions 
of the Mott phase boundaries. Consider therefore the 
particle-hole symmetric model (j2.4p with, for concrete- 
ness, Jij = Jo{l + AJij) > 0, nonzero on nearest neighbor 



C. Correlations and excitations in the random 
bond model 

Consider a superfluid droplet with (spatial) volume V, 
which will occur roughly with density e~^°^, for some 
constant po- The behavior of x cxd cylinders of magnet 
has been discussed in detail by Fisher and Privman^ 
who were concerned with finite size scaling theory of fer- 
romagnets with a continuous 0(ri) symmetry below their 
bulk critical points. Their main result (which will be 
rederived in a more general context below) was that the 
correlation length, ^||, along the cylinder is governed by 
the bulk helicity modulus along the same direction: 



^11 



2T(T)y 
{n~l)kBT' 



(3.15) 



In our case, ksT = 1, 71 = 2 and T(T) = T^(Jo). The 
correlation function, Go(t), along the cylinder then varies 
as 



Go(r) 



(3.16) 
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There is some ambiguity in what we should take for V 
and Tt{J) in (|3.15p : the droplets are neither perfectly 
spherical, nor is Jy uniform throughout the droplet. 
Thus V should be some effective volume, while Tt-(Jo) 
should be the bulk temporal helicity modulus associated 
with some effective uniform Jq > J^, say roughly the av- 
erage of Jij over the droplet. None of these ambiguities 
change the order of magnitude estimates made below. 

1. Stretched exponential correlations in the symmetric glass 

The full temporal correlation function, G(t), is ob- 
tained by averaging Go{t) over all droplets (considered 
independent in the present picture). We therefore esti- 
mate 

G(r) ^ JdV J dTrp{V,Tr)Go(T), (3.17) 

where p{V, T^) is the probability density for droplets of 
volume V and bulk helicity modulus T,-, 

p{V,Tr) ^ e-^/^°^^-K (3.18) 

The coefficient Vq{Tt-), which we interpret as the "typi- 
cal" droplet size for a given T^, will depend on the de- 
tailed shape of the tail of the probability distribution for 
Jij > J°. Using (|3.15[) . for large r we may perform the 
integral over V using the saddle point method. The in- 
tegration will be dominated by V near the solution of 
^(WK) + T/2TrV) = 0. This yields 

G{t)^ j dT,e-V2-/T.Vo(TO^ T^oo. (3.19) 

The coefficient Tt-Vo(Tt-) will have a minimum at some 
value, T,-(Jo), corresponding to the most probable large 
droplets, and this will govern the asymptotic behavior of 
the integral (|3.19p to yield finally, 

G(T)^e-V-/-o(^o), ro(Jo) = ^T,l^o(TO. (3.20) 

The droplets therefore yield a stretched exponential be- 
havior, to be contrasted with the purely exponential be- 
havior in the Mott phase,^ From ([X^ we may de- 
rive the quantum mechanical single-particle density of 
states, -J^^ pi{e), defined as the inverse Laplace transform 
of G(t): 

G{t) = \ depi(e)e"'l^l. (3.21) 

It is easy to see that exponential decay in G(r) requires 
a gap in p\[e), 

pi(e) = 0, e<e, ^ G(t) ^ e-^=l"l, (3.22) 

while slower than exponential decay permits p\(t) > 
for aU e > 0. The form ([3:20ll yields 



a "soft gap." The non-exponential form of (|3.20[) and 
the gap free form of p.23p are known as Griffiths 
singularities,— and define the RRG phase for Jo within 
some interval above Jq < J° — Aj < Jc- 

2. Lack of a direct Mott-super fluid transition 

The fact that there cannot be a direct Mott-SF tran- 
sition (i.e., Jc = Jq) follows from the fact that the corre- 
lation lengths in both the superfluid droplets and in the 
"background" Mott phase (defined, for example, as the 
region between droplets in which all Jij are some speci- 
fied finite distance below J°) are finite. 

The superfluid transition can occur only if the droplets 
grow to be large enough, and/or close enough together, 
that they begin to coalesce. Thus, we have treated the 
droplets as independent, but there will actually be expo- 
nentially small interactions ~ Q-diJo)/^(Jo) between them 
due to the finite background correlation length, where 
rf(Jo) is the typical droplet separation [which diverges 
exponentially as Jq Jc/i^ + ^+)] and C(Jo) is the ef- 
fective correlation length in the surrounding Mott phase 
(which remains finite in this same limit). These cou- 
plings will increase the correlation length slightly, but 
only when Jq is sufficiently large [a finite distance above 
J°/(l -I- S+)], and do (Jo) sufficiently small, will the su- 
perfluid droplets effectively overlap, and the correlation 
length diverge. This defines Jc, which, as noted earlier, 
could lie below J° for certain disorder distributions. 

3. Correlations, excitations and compressibility at small 
nonzero fi: Bose glass onset 

Now consider the compressibility. Its computation re- 
quires the addition of a small uniform chemical potential, 
fi. As alluded to earlier, we expect pi (e) to be finite at 
e in the presence of ^, signifying a Bose glass, and 
implying power law behavior for G(r)i^ 

G(r) «pi(0;Ai)/r, r ^ 00, (3.24) 

though we shall see that pi{0; p) is exponentially small 
in i. 

Such behavior lies far outside any classical intuition. 
To see how it comes about we must generalize the ideas 
of Ref. [33| to this case. Fortunately this is relatively 
straightforward: a compact statement of the Fisher- 
Privman result is that long time correlations along V xoo 
cylinders (for n = 2) are governed by an effective one di- 
mensional classical action 

^Off = -^VTr dT[dr<P{T)]\ (3.25) 

where 4'{t) is a coarse-grained phase. This immediately 
yields 

Go(r) EE (e'[*(-)-*(0)l) = e-"^mr)-m?) ^ (3.26) 
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which, upon using 



1 



1 



duj 



, tj2 2TrV' 

(3.27) 

yields ((31^ and ((3A6)) . 

Now we must generahze (I3.25P to finite /i. This is 
accomphshed using p.6p : effective long wavelength, long 
time "hydrodynamic" fluctuations in the phase cf) are gov- 
erned by precisely the same elastic moduli that govern 
equilibrium twists in the phase. Thus in p.6p one sim- 
ply replaces by 9^0 and integrates over all space. If, 
as in the present case, the twists in different directions, 
a, linearly superimpose [this may be checked directly 
from (13. 7p . where now one defines = 6^'^°*°^°/^°?/;, 
with ij) obeying 0Q,-boundary conditions simultaneously 
in each direction], one simply sums over all directions a 
to obtain the final result: 



5e 



D 

E 



d'^xdr 



(3.28) 

In the case where the interactions are spatially isotropic 
one has pa = and Tq = T for a = 1, . . . , d. With the 
identifications (j3.13p for a = t we have 



d'^xdr 



-ipdrcf) + ^K(5.r(/ 



(3.29) 

The validity of this hydrodynamic form in the presence 
of disorder relies on a hidden assumption that no new, 
unforseen low energy excitations develop, e.g., in the am- 
plitude, rather than just the phase, of the order parame- 
ter. This appears unlikely, and there are concrete propos- 
als for such excitations (but see Ref. [l3 for some discus- 
sion on this point in the context of interpreting quantum 
Monte Carlo data). 

For X oo cylindrical geometries, the effective one di- 
mensional result, (|3.15p and p.l6p . is obtained by as- 
suming that for each r, (/>(x, r) is essentially constant in 
space, and hence that only the temporal fiuctuations are 
important. More formally, the finiteness of V implies an 
energy gap in the spatial spin-wave spectrum, between 
uniform 0(x) and the next excited state in which (j) twists 
by 2tt from one side of the system to the other, of order 
The temporal spectrum has no such gap (the 
frequency, uj, in (|3.27p is continuous), and therefore the 
asymptotic long time, large distance behavior may be ob- 
tained by assuming (/)(x, r) — (pir) only. The |V(/)p term 
in p.28p may be treated as an additive constant that 
drops out of any temporal average, and we obtain the 
proper generalization of (I3.24p : 



Let us now study the consequences of (I3.30p . First, 
when pV is an integer (i.e., the density in the bulk is com- 
mensurate with the droplet volume, V) the 27r-periodic 
boundary conditions on (j) imply that the p term sim- 

ply drops out of the statistical factor, e'^^^ff . This im- 
plies that only the fractional part, pV mod 1, matters 
in (j3.30p . One must be careful to distinguish p and k from 
the actual density and compressibility of the droplet of 
volume V . The values of p and k, are appropriate to a 
bulk superfiuid system with some effective J > J^. The 
bulk compressibility, kq(J), of such a system is finite and 
nonzero. When p = the density is p = (or, more 
generally, some integer), so for small p, 

p = Ko{j)p + O{p^) 

K = Ko{j) + O{p). (3.31) 

The actual values of p and k in the droplet are now 
computed from (|3.30p as follows: the free energy density 
is given by 



(3.32) 



k{Po) 



in which /o is the bulk free energy density corresponding 
to the input parameters, k and p. Thus, for example, at 
a given value of the chemical potential, p = po, we have 

(3.33) 

and hence, correct to quadratic order in p — po, we may 
take 

foip) - foM - P°(/i - Mo) - ^«°(/i - Pof. (3.34) 

The effective action must also be correct to quadratic 
order, therefore for the purposes of computing the full 
free energy, consistency requires that in we take k = 
Hp and p = p^+K^{p—pQ). For the purposes of computing 
derivatives with respect to p, the only /^-dependence in 
the fluctuation part of the free energy is now in p. We 
obtain 



OO 

E 

_m— — o 



E 



tr* 



(3.35) 



where tr"^ means that we impose the temporal boundary 
condition (/)(/3) — (t){0) + 27rm. Now define (?i(r) — (pir) — 
27rTOr//3, so that — '?i(0), to obtain 



E 



e ^ e 



'eff 



r0 

-V I dr 
Jo 



]^Hi{dr(j)f' - ipdrCj) 



(3.30) 



All the effects of particle-hole asymmetry are in the p 
term. 



X tr" 



= /o + /oo['*1-^ln 



E 

./ — — OO 



(3.36) 
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where we have used (see App. |B|) 



E 



m— — oo 



with K = /3/47r2K0\/, and 
/oo[«°] = In 



E 



-27r^_R'(a;-/)^ 



(3.37) 



(3.38) 



is independent of fi — fiQ. In the hmit P —^ oo only the 
term with minimal {x — Z)^, i.e. — ^ < a; — Z < 5, con- 
tributes (at the boundaries, two neighboring terms are 
degenerate). Let Io{k'^,Ijl) be this minimizing value of I. 
We then obtain finally, 

/ = ./o + ./oo + ^(pl^-W', (3.39) 
and the actual density in the droplet is 

There are exactly Iq particles in the droplet for the in- 
terval of /i such that \p{iJ,)V — lo\ < ^, and we have 
established the desired result that the droplet is incom- 
pressible on this same interval. 

Consider next the temporal correlation function, given 

by 



tr-^ 


e ■=« e 


»[0(r)- 


-Hr')] 


tr"^ 







(3.41) 



l^m=-oo ^ 


o(0) 

e'^=ff e 


i[0(r)- 


-0(^')]" 




o(0)' 





Defining the same periodic field, (/'(t), we obtain 
G(0)(T-r') = (e^['^M-'^(^')l)„(o, 



EDO 
m— - 



^m— — 00 

Hr-T'|/2K''yg-(r-r')(pV modl)/«:OV^ /3 ^ OO, 



(3.42) 



where we have used (I3.37p . Once again, in the limit 
/3 ^ 00 only the term with — ^ < pV — I = pV 
mod 1 < \, contributes. One sees now that Gp{T) de- 
cays exponentially for both r ^ ±00, but at different 
rates: 



G(")(t) = e-(i±7)|r|/2KV^ r^±oo 
-1< 7 = 2(py mod 1) < 1. 



(3.43) 



This exponential decay signifies an energy gap, propor- 
1 

W 



tional to tw, for adding a particle, and is equivalent to 



the incompressibility result above. However, for large V 
this gap is very small, and one need only increase p (and 
hence p) by a small amount to add a single particle to 
the droplet. 

For given p the number of particles in the droplet will 
be I = [pV], the greatest integer less than or equal to pV . 
Since there exist arbitrarily large droplets, an arbitrar- 
ily small change in p will add particles to the system in 
precisely those droplets with volume V > ^ k, Fo- 
cusing on p near zero (where 
the total density as 



1 

K.°^J^■ 

Ko), we may estimate 



Ptot 



dVdKop{V, Ko 
Kop 



[Ko^J■V] 
V 



Kope 



"Of 

-l/SoMVo(fto) 



(3.44) 



where Rq is defined analogously to T.^ in (|3.20|1 . In the 
derivation of this formula we have assumed that p > 0, 
but the result is valid also for /i < if is replaced by 
in the exponent (only). The total compressibility may be 
estimated as 



Utot 



1 



Vo\^^\ 



,-l/KoVb(Ko)lM 



(3.45) 



which also vanishes exponentially as \p\ —> 0. 

Finally we may use the above results to estimate the 
total temporal correlation function and to exhibit the 
finite density of states, p.24|) . at e = 0. Once again, the 

total correlation function, Gpij) is the average of Gp°^ (r) 
over all droplets: 



Gp(T) = / dVdKop{V,Ko)G'^;^\T;Ko,V). 



(3.46) 



For large r and small p > 0, only large volumes con- 
tribute to the integral. It is clear from p.43p that Gp^-* (r) 
decays most slowly when pV is close to half-integer, and 
those droplets with such "resonant" values of V will con- 
tribute the leading large r-dependence. The smallest res- 
onant volume (into which a single particle will be added) 
is precisely V = and contributions from higher order 



resonances, V = 



., will be exponentially smaller 



in -. Thus 



Gp{T) - / dT/dKop(l^,Ko)e-l"'/'"°^e-^"/2'^°^ 

dx 



-l/2Ko|M|Vb(Ko) 



'^olMl 



(3.47) 



where x — \pV — i| and 5 < ^ is a cutoff and we have 
replaced V by its smallest resonant value, ^ « ji^' 
cry where except in 7 = 2{pV mod 1). The integration 
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FIG. 5: (Color online) Schematic illustration of the droplet model of the Bose glass phase for the random site energy model. Left: 
Exiting the Mott lobe at fixed Jo < Ja' below its tip with increasing > /i+(Jo) — A [or decreasing < /i_(Jo, A) = /i_(Jo)+A]. 
For arbitrarily small e = /i — (Jo, A) there will be an (exponentially small) density of droplets with volume large enough that 
the energy gap to adding a particle is smaller than e. Right: Exiting the Mott lobe from its tip at fixed = with increasing 
Jo > Jo'', in the statistically particle-hole symmetric model. For arbitrarily small A J = Jo — Jo^(A) there will similarly be 
arbitrarily large droplets whose local chemical potential magnitude is above the energy gap required to add a particle (previous 
droplets from the left illustration), or remove a particle (additional droplets). 



is now trivial, and we obtain 



Gp(T) 



■ e 1 



S\ut\ 



^-1/2ko|/j|V'o(ko) 



(3.48) 



This reproduces the ^ behavior (I3.24p . predicted for the 
Bose glass phase with 



Pi(e-O) 



-l/2Ko|Ai|Vo(Ko) 



(3.49) 



Note that the power law prefactors (in fx) of the exponen- 
tial must not be taken seriously because we have made a 
very crude estimate for the probability function p{V^ kq)- 
Recall that kq is the "most probable" compressibility for 
large droplets. 

One direct consequence of the slow power law decay 
of temporal correlations, (|3.24|) or p.48p . is a divergent 
superfluid susceptibility)^ 



Xs 



j (fx j dTG{x,T) 

Voifia) j drGpir)^^. (3.50) 



This provides another signature, in addition to the finite 
compressibility, distinguishing the Bose glass from the 
Mott phases. 

To summarize, we have seen that for the particle- 
hole symmetric model the correlation function, G(r), has 
stretched exponential behavior coming from large rare 
regions in which J > J^. This is known as a Griffiths 
singularity^- , and this kind of effect is ubiquitous in ran- 
dom systems. Since G(r) still decays faster than any 



power law, the effects of these singularities are obviously 
physically rather subtle. In contrast, when fi ^ the 
model no longer has a classical interpretation, and the 
behavior is far more singular: for given fi, finite droplets 
of size V ^ ^Kol/xj give rise to power law decay of G'p(r) — 
no longer do the singularities occur only in the limit 
^ oo. Quantum mechanically, we understand this 
as being a consequence of the existence of arbitrarily low 
energy single particle excitations, arising from superfluid 
droplets with very small energy gaps for the addition of 
an extra particle. It is interesting to see this derived ex- 
plicitly from the interference terms in the Lagrangian [see 



D. Droplets in the random site energy model 

1. MI-BG phase boundary 

Consider now the random site energy model, with uni- 
form (nonrandom) hopping Jij. Above, we studied the 
crossover between the RRG and BG phases with applica- 
tion of a small uniform with Jq beyond the tip of the 
Mott lobe. Here we begin by considering Jq < Jo^(A) 
below the tip of the Mott lobe (to be computed below), 
and consider, for specificity, the transition to the BG 
phase with increasing ^ (identical arguments, with parti- 
cles replaced by holes, go through for the transition with 
decreasing p). 

For ji < /i+(Jo) — A, all local chemical potentials fi — Si 
lie below the Mott gap /^+(Jo), and no extra particles 
can enter the system (if a uniform chemical potential /i 
lies below the Mott gap, then reducing it on some sites 
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Droplet model schematic 
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FIG. 6: (Color online) Schematic illustration of spectrum of 
droplet excitations. Droplets are assumed to be well sep- 
arated, and independent (see Fig. [S]). The horizontal axis 
is some convenient droplet index. The vertical axis shows 
the critical chemical potential levels (in some arbitrary units) 
at which a new particle is added for each droplet. The 
lower bound of the excitation spectrum lies at the Mott 
phase boundary fi+{Jo) — A (solid horizontal line). For any 
e = /i— /i+(Jo)-|-A > (dashed horizontal line), there will be a 
finite (but exponentially small) density of large droplets with 
local chemical potential lying sufficiently far above /^+(Jo) 
that one or more extra particles are added. The total number 
of particles added is given by counting the number of "occu- 
pied" levels below /i. For a given droplet, while the lowest 
excitation depends primarily on the local chemical potential, 
subsequent particles are added in sequence, with gaps scaling 
as IjKj^V where V is the droplet volume, and K+(Jo,/i+) is 
the bulk compressibility of the pure (superfluid) phase just 
above Mott phase boundary. In an infinite system there are 
infinitely many droplets, implying a continuous distribution 
of excitation energies, and the bulk density will increase con- 
tinuously with increasing e > 0. 

cannot reduce the gap'^'*). Thus, /i+(Jo) — A is a lower 
bound for the Mott phase boundary. 

Another large rare region argument now shows that it 
is also an upper bound. Let e = — /i+(Jo) -I- A > 
be positive, but arbitrarily small. Then there will be a 
finite density ~ e"^/^'^^-* of arbitrarily large, but very 
rare regions of volume V, in which all of the site energies 
are larger than, say, A — e/2 — see the left panel of Fig. [5l 
The scale Vb will diverge as e ^ in a fashion determined 
by the shape of the tail of the distribution. The region 
will then look like a finite size droplet with local chemical 
potential larger than /i+(Jo) by at least e/2 [additional 
near-uniformity constraints may be imposed on the site 
energies if necessary to make the droplet as homogeneous 
as desired, but this will not be critical to the argument as 
it entails only an adjustment of the definition of Vb(e)]. 

Let K_|_(Jo) be the (finite) compressibility of a homo- 
geneous system just above the Mott lobe. From the Bo- 
goluibov theory of the dilute Bose gas of quasiparticle 
excitations, one may estimate k+ ~ m/?i^^(Jo)''~^ (re- 
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FIG. 7: (Color online) Mott lobe boundaries as a function of 
disorder. The phase boundary is determined by that of the 
pure system and the half-width A of the disorder distribution: 
/i±(Jo, A) — /x±(Jo, 0) =F A. The tip of the Mott lobe, where 
a statistical particle-hole symmetry obtains, occurs at Jo = 
Ji'iA) defined by ^±(Jo) = A. 



placed by a logarithmic form in d = 2^^), where £,{Jo) 
estimates the diameter of the quasiparticles which de- 
termines the s-wave scattering length.— Then the added 
number of particles may be estimated as I = [eK+V/2] 
[compare p.43p and below]. Thus, if > 2/K+e the 
droplet will have at least one extra particle, and addi- 
tional particles are added each with energy gap 2/k+V. 
Note that one should have as well V :s> ^( Jo)'', the quasi- 
particle volume, for this dilute Bose gas argument to 
make sense. The droplet excitation spectrum following 
from these arguments is illustrated in Fig. |6l 

A calculation identical in form to ()3.44|) and ()3.45p 
(with /i replaced by e/2) can now be used to show that the 
bulk compressibility is finite, though exponentially small 
in 1/e. This finally demonstrates that fi = /x+(Jo) — A 
is also a lower bound on the Bose glass phase boundary, 
and hence is, in fact, the phase boundary. The resulting 
shrinking of the pure system Mott lobe is illustrated in 
Fig. [7] for various values of S = A/Uq- Note that the 
pure system fJ,±{Jo) ^ [Jo — critical singularity^ is 

replaced by a slope discontinuity at the tip of the Mott 
lobe, Jo = J^^ [defined by n±iJi^) = A]. 



2. Statistical particle-hole symmetry 

Since we assume that the random site energies have a 
symmetric distribution, the line fi = 0, passing through 
the tip of the Mott lobe, has a statistical particle-hole 
symmetry (see the discussion in Sec. Does this 

change the nature of the Bose-glass phase along this line, 

Jo > 4'^ 

It is clear that the answer must be no: for any 
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Jo > Jq^ , hence Sn = A — |/x±(Jo)| > 0, we will be 
able to find an interpenetrating, but independent, dis- 
tribution of droplets of arbitrarily large size in which 
all Ei < — /i_|_(Jo) — or all > — ^_(Jo) + Sfj./2. 

The former will have additional particles, while the lat- 
ter will have a particle deficit (additional holes) — right 
panel of Fig. [5] Precisely at = 0, symmetry implies 
that the overall density is still fixed at the Mott value. 
The previous analysis, generating finite compressibility 
and nonzero energy density of excitation states, then goes 
through precisely as before, but now with e replaced by 
6^, and a different volume scale, Vq{k±, now depend- 
ing on 

Thus, at least at this qualitative level, statistical 
particlc-hole symmetry differs from generic particle-hole 
asymmetry only in that there are now two sets of droplets 
(particle droplets and hole droplets) contributing inde- 
pendently to the excitation spectrum. It seems very un- 
likely then that the nature of the superfluid transition 
would be any different either. We shall address this issue 
further in Sec. IIVI 



3. Lack of a direct Mott- superfluid transition 

Using arguments similar to that presented in Sec. 
nil C 2[ we can rule out a direct MI-SF transition in this 
model as well. The excited droplets, containing a di- 
lute superfluid of excitations, have been treated as inde- 
pendent. There will, however, be exponentially decaying 
interactions ^ Q-d/HJo) between them, where d is the 
typical droplet separation, and ^(Jq) is the correlation 
length in the background insulating Mott phase. Even 
for very large d (exponentially large in 1 /e) , it is not im- 
mediately obvious that there might be some (exponen- 
tially) small hopping of quasiparticles between droplets, 
generating bulk superfluid coherence. However, the exci- 
tation spectrum of each individual droplet is discrete, and 
the usual Anderson localization arguments imply that for 
small e the competition between hopping distance and en- 
ergy level matching implies that the excited states remain 
strongly localized to the neighborhood of their droplets. 
Stated shghtly differently, for e not too large, the excited 
particles see a residual random potential whose effective 
low lying single particle states must be localized. Identi- 
cal arguments, applied to the two sets of droplets, imply 
that there can be no direct transition through the statis- 
tical particle-hole symmetric tip of the Mott lobe either. 

In apparent violation of these, essentially rigorous, ar- 
guments, there have been some recent Monte Carlo sim- 
ulations claiming to see a direct MI-SF transition, over 
a flnite segment of the Mott lobe around /U = 0, for suf- 
ficiently weak disorder and even claiming evidence for 
new multicritical behavior at the endpoints of this seg- 
ment. However, it is well known that rare region effects 
are invisible to Monte Carlo simulations, ^'^^ which are, 
by necessity, limited to finite volumes with at most a few 
thousand sites. 



Note, in addition, that since the boundary of the Mott 
phase is known exactly, a direct transition would imply 
that the SF phase moves in to take over all of what used 
to be the pure system Mott phase. Although disorder 
can indeed increase the stability of the supcrfiuid over the 
MI,^^ that it would eat up the putative intervening BG 
phase entirely is disproven by the rare region arguments. 

At the risk of belaboring the point, notice also how ex- 
traordinarily sensitive to the type of disorder the direct 
MI-SF transition would have to be if it existed. Con- 
sider a model in which the usual top-hat disorder model 
of sufficiently small half-width A is augmented by an ad- 
ditional Gaussian disorder with the same width A, but 
with a miniscule relative amplitude — say 10~^. Since the 
onsite potentials are now unbounded (though one would 
have to survey billions of sites to know it) the MI phase 
no longer exists. However, by any conceivable measure, 
the total disorder is still small, but a direct MI-SF tran- 
sition would now require a SF phase all the way down to 
Jo = 0. In other words, a one-in-a-billion change must 
cause the phase boundary to move an infinite distance 
(on the natural scale of I/Jq). It is obvious, however, 
that the simulations would see no change at all with the 
infinitesimal added Gaussian. 

Regarding the apparent multicritical scaling, it is very 
easy to find apparent scaling of data over the limited 
ranges of system sizes that are available, if one has a 
sufficient number of free parameters available to fit (the 
transition point, and the exponents v and z in this case). 
A more likely explanation for the apparent scaling is a 
combination of finite size effects, and a crossover from 
the pure system critical behavior to dirty boson critical 
behavior at weak disorder. The latter implies a crossover 
scaling variable of the form A/\Jq — J^\'^, where (j) is 
a crossover exponent quantifying the instability of the 
pure MI-SF transition to small disorder. At small ^, this 
crossover will also mix with the pure system particle-hole 
symmetry-breaking scaling variable fJ./\Jo — Jc\'^° ■ For 
small A and limited system sizes, these scaling variables 
will saturate before the asymptotic BG-SF dirty boson 
criticality can become visible very close to Jo,c(m)- The 
corresponding multi-crossover scaling form, which will be 
discussed in more detail in Sec. lIVBl below, could easily 
mimic multicriticality. 



IV. PARTICLE-HOLE SYMMETRY AND 
SCALING NEAR CRITICALITY 

In order to discuss scaling it is convenient (but by no 
means necessary) to use the Lagrangian, (|2.10p . fur- 
ther simplified by taking the continuum limit and drop- 
ping all unnecessary dimensionful coefficients. To begin 
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we take Cr = 1 only, and consider the Lagrangian, 



+ ir(x)|^p + iu|^|4 



(4.1) 



equivalent to £5 in Table U with g{x.) — go + Sg{x.) and 
r(x) = ro + Sr{pc). The phase transition occurs when the 
control parameter ro becomes sufBciently negative. As 
described earlier, when g = the particle-hole symmetric 
problem is recovered. If the IStiAP term is dropped and 
we take g = I, we obtain the closest approximation to 
the boson coherent state Lagrangian, (12.11) . which was 
the starting point for the work in Ref. We shall see 
that the term, which was ignored in Ref. is 

actually crucial for a correct understanding of the critical 
behavior. When g = one obtains precisely the random 
rod model studied in Ref. [l^. 

As part of our scaling discussion, we will revisit pre- 
vious arguments^ for the scaling relation z = d for the 
dynamical exponent for the dirty boson problem. Re- 
cent quantum Monte Carlo results contradict this rela- 
tion, finding z = 1.40±0.02 in d = 2^ We will show that 
all of the previous arguments, when looked at more care- 
fully, in fact place no constraint on the exponent zi^ 
Lacking deeper arguments, it would appear that z re- 
mains an independent exponent, undetermined by any 
scaling relatioui^^ 



A. Scaling of superfluid density and compressibility 



Recall that, as described in Sec. IIII A[ the superfluid 
density, ps^ and compressibility, k, measure the system 
response to twists, spatial and temporal respectively, in 
the superfluid order parameter. As in (j3.3p . we introduce 



(4.2) 



V;(x,t) = e^*('^«-''+'^''^V(x,r), 



and impose periodic boundary conditions on ?/'(x, t), 
identifying LOa = 60/ (3 and ko = {9i/ Li, . . . ,9d/ Ld). 
Analogously to p.4p . one obtains 



(4.3) 



with 



5CM = - J d'x I dT{^ikl+u;lM' (4.4) 
- «woV'*[c't - .9(x)]?/' + iko • ?/'*V?/'|. 

The expansion of the free energy p.Sp in powers of ko 
and LOo takes the form 

5f{ko,uJo) = -ipujo + ]^{ni^l + ^kl) + 0{u}l,ujokl), (4.5) 



In addition to the result (|3.14p for the helicity modulus, 
one may then identify 

K EE Tr=[{\m],,+ J d'x J dT 

X [ir {dr - 5) V] (x, r) ir (dr ~ g)i'] (0, 0)),],, 

p ^ Pr = -[{ridr~gmi,. (4.6) 

The spatial isotropy of (|4.ip implies that = T is inde- 
pendent of the direction a. The subscript "c" on the right 
hand side of the second equation indicates a cummulant 
average, i.e., that the product of the averages, namely p^, 
should be subtracted from the integrand. We shall ulti- 
mately require these expressions only when g = 0, where 
p = as well. 



1. Random rod critical point scaling 

Let us first consider the scaling of ps and k for the clas- 
sical random rod problem, g{x) = 0. Note that ko-V'*V'0 
and ujo4'*9t'<P are symmetry breaking perturbations: the 
first breaks the x ^ — x spatial inversion symmetry, and 
the second breaks the time inversion symmetry t —t 
symmetries of random rod Lagrangian. The latter corre- 
sponds precisely to particle-hole symmetry. 

Critical universality classes are, by definition, insensi- 
tive to most changes in the detailed parameters of the 
Hamiltonian, but symmetry breaking perturbations are 
often an exception, leading to changes in the asymptotic 
critical behavior. One therefore expects fcp, loo to be rele- 
vant perturbations to the random rod problem, entering 
the thermodynamics through scaling combinations that 
diverge as the critical point is approached. Since fco is 
an inverse length and uio is an inverse time, one expects 
them to be scaled by the corresponding divergent corre- 
lation length and time, respectively. This motivates the 
following form for the singular part of the free energy: 



/,(fco,c^o)« A|5|2~«ci>(fcoe,u;oer), 



(4.7) 



where ^ « ^o|<^| and ~ ^t,o|'^|~''^" are the cor- 
relation lengths in the spatial and temporal directions, 
respectively, A is a nonuniversal amplitude, and the dy- 
namical exponent is defined by zq = Vro/vo- The sub- 
script on the exponents indicate that they are those 
appropriate to the classical random rod problem, and 
the generic auxiliary parameter, 5 is ro — ro,c in (|4.ip . 
but more generally is any parameter such as chemical 
potential, pressure, strength of disorder, film thickness, 
or magnetic field, which moves the system through the 
phase transition at T = 0, defined to occur at (5 = 0. We 
assume that (5 > corresponds to the disordered phase 
and (5 < to the ordered (superfluid or superconduct- 
ing) phasei^l The two scaling arguments fcoC and wo^t 
indeed diverge as ^ — *■ for arbitrarily small ko and loo^ 
consistent with their expected relevance. 

Although (|4.8p motivates the correct result, the fact 
that LJo and fco are infinitesimal in the zero temperature 
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thermodynamic limit, (3, La — > oo, means that a httle 
more care is required to construct a rigorous scahng for- 
mulation. More properly, the boundary condition depen- 
dence appears in a finite size scaling ansatz for the fi'ee 
energy;^ 

5f « (3-'L-''^l[5{L/iaY'''-,5{P/^r.»f'''-°]. (4.8) 

The existence of a nonzero stiffness (in the ordered 
phase), i.e., via (|3.6p . a leading finite-size correction of 
order or now requires that the scaling function 
obey ^l{x,y) « a;rf''oyZoi'o($e^-2i.o _^$ey-2zoi^o) foj. j^rge 

(and 5 < 0), yielding 

T ~ e^''C(l($?/2^202)^.o 

^ ~ £.^%fi{^ll2ii^e^)5^^<^ (4.9) 

with the Josephson scaling relations,— 

vo {d + za - 2)iyo ^ 2 - ao - 2vo 
VtO ^ [d- zo)iyQ = 2 - ao - 2zqVo, (4.10) 

and requiring in addition that ^12'^^^- We emphasize 
that the crucial assumption is that the leading boundary 
condition dependence is all in the singular, i.e., finite 
size scaling part, of the free energy. We can make this 
assumption only because ko and ujq introduce relevant 
perturbations which fundamentally alter the symmetry 
of the Lagrangian. Further support for this assumption 
is that we expect all stiffnesses to vanish identically in 
the disordered phase of the classical model: thus T and 
K can have no analytic contributions at all. 

2. Dirty boson critical point scaling: revisiting z — d 

Now, we can try to extend the above arguments for 
nonzero following Ref.@. Thus, one may posit identical 
forms (14. 7p or (|4.8p for the singular part of the free energy, 
but with exponents and scaling functions appropriate to 
the dirty boson critical point. One obtains then 

T - |(5|'', v = 2-a-2v={d + z-2)v 

\5\''^, Vr = 2-a~2zv^{d~z)v. (4.11) 

For 5^0, both the Bose glass and superfluid phases are 
compressible, so it is expected that the compressibility 
remains finite right through the transition. This leads to 
the prediction z — d. 

However, although the result for T is believed to be 
correct, there are a number of questionable assumptions 
underlying the argument above for k, as we shall now 
discuss. For g ^ Q, the density of the Bose glass phase 
varies smoothly with (the control parameter analogous 
to the chemical potential ji) . A temporal twist only per- 
turbs slightly the term, gip*drtp, that is already present in 
the Lagrangian, and is therefore not expected to produce 
a new relevant perturbation. Thus, the scaling variable 
cannot be wqCt- Rather loq produces only an infinitesimal 



shift go ^ go ~ i^o^ which, through the above analytic- 
ity argument (including also the possibility of complex 
shifts), alters the free energy in a completely predictable 
fashion unrelated to scaling. By way of contrast, spa- 
tial twists still represent a relevant perturbation, and we 
therefore predict a scaling form at tjQ = 0: 

5f - r 'i-'<I''[5(i/eo)'/^ (5(/3/C..o)'/''1, (4.12) 

with (^{x^y) « ^-^x^d.-2)uyzu }g^j.gg 2; yielding T « 
($^/27r2e'2)J", with V ^ {d + z-2)v = 2 - a-2v &s 
before. We expect only small subleading corrections in 
y^'^: since no temporal twist has been imposed, there can 
be no 0(/3~^,/3~^) corrections. 

Now, if we include a finite wq, the basic change in (|4.12p 
is that go 9o — i^o everywhere, and in addition one 
must include changes arising from boundary condition 
dependence of the analytic part of the free energy. One 
obtains 

5f = /3-iL-''<f^[<5,(L/Co)^/^<59(/3/Cr,o)'/''1 

+ /°(ro,.9o-*^o)-/°(r-o,.go), (4.13) 

where /° is the free energy density in the absence of all 
twists (i.e., under fully periodic boundary conditions), 
including both analytic and singular parts, and Sg = ro — 
ro,c{go — i'-^o) ~ S + iuJoTQ c(ffo) is the perturbed deviation 
from the critical line ^^0(30)- Most importantly, <i>^ is the 
same function as that in (|4.12p . and therefore produces 
only small corrections in (/3/Ct)~^- The scaling function 
itself therefore contributes only vanishing corrections to 
K in the limit /3 — > cxd. 

All contributions to k are therefore contained in f^, 
and arise from (a) its analytic part, (b) the (linear) ujo 
dependence of 5e in its singular parti^S Regarding (b), the 
leading singular part of /° is of the form /s^^g ~ 
The go-dependence of ro^c couples derivatives with re- 
spect to go to those with respect to ^, and one obtains 
therefore contributions psing ~ |<5e|^~" to the density and 
Ksing ^ I'^el " to the compressibility. The hyperscaling 
relation a = 2 — {d + z)v implies that a < under the 
rather weak condition that v > 2/{d + z) (which by all 
evidence to date appears to be strongly satisfied)^) and 
the singular parts of p and k are expected to vanish at 
criticality. Put slightly differently, given that n is already 
nonzero in the Bose glass phase, signifying the existence 
of long range temporal correlations on both sides of the 
transition, it would not be surprising to find that the 
critical singularity, signifying the adjustment of the den- 
sity of states pi(e = 0) as the density of large droplets 
increases, leads to only small corrections to n. 

Regarding (a), the analytic part of the free energy has 
an expansion 

fa{ro,9o) = -Pc(ro)[ffo-go,c(ro)] (4.14) 

- ■^i^c{ro)[go - goA'^o)? + 

about the transition line go.cifo)i where Pc{fo) and Kci^o) 
are now recognized as the finite values of p and k at the 
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transition. These finite values are predicted irrespective 
of the value of the exponent z, now seen as unconstrained 
by the present arguments. Note further that since the 
leading cjo-dependence is linear, if the leading bound- 
ary condition dependence were indeed through the scal- 
ing combination luq^t (or, more properly, the finite size 
scaling variable 5/ [3^^^"), the leading contribution to the 
density must take the form p ^ \5\'^'^ ^ contradicting the 
fact that the density must be finite at the transition. 

To end this discussion, it is worth emphasizing why the 
arguments above fail for the classical random rod model, 
f;(x) = 0. The key point is that ro,c(5o) is singular at the 
special value go = 0, and 5g is no longer an analytic func- 
tion of Wo- Essentially, the special symmetry at 50 = 
implies that 5 and /i are "orthogonal" thermodynamic 
coordinates and the derivatives with respect to /i that 
define k = do not mix with derivatives with respect 
to 5. We have already seen, therefore, that k = in the 
disordered phase. We therefore expect k, to rise contin- 
uously from zero for (5 < 0, with the exponent Vr > 0. 
This implies that z < d in this case (equality is still per- 
mitted and would imply a discontinuity in ft at (5 = 0, 
which indeed is the case in d = 1 — see Sec. |V]). Note 
that for homogeneous classical disorder, where the coef- 
ficient r in (|4.ip depends on both x and r, we will have 
isotropic scaling, z = 1. The rod disorder should increase 
z. 



3. Scaling of correlations 

Let us now consider the two-point correlation function, 

G(x,r) = [(V^(x,r)^(O,0)),],,, (4.15) 

whose Fourier transform is normally assumed to scale in 
the form. 



If one naively matches (|4.16p and (|4.17p . one concludes 
that g{x, y) « {gix^ + gilj^) for small x, y and hence that 



9ia 



(4.19) 



Using the well known scaling relation a -f 2/3-1-7 — 2, one 
recovers (|4.11l) . 

However, (I4.19P is really just a disguised version of the 
free energy argument. Thus, ScS assumes that the ener- 
getics of global phase twists also describes slowly varying 
local phase twists. Thus, locally we replace ko by V(/) 
and Wo by drfj), then integrate over space-time. Once 
again, (|4.16p - (|4.19p are expected to be valid for the ran- 
dom rod problem. However, for the dirty boson problem, 
if K arises from the nonscaling part of the free energy, it 
is unlikely that it can now arise from the scaling part 
of the two-point function. Thus, in place of (|4.16p . we 
propose instead the scaling form 



G(k,w) 



D\d\^-'-TiH,LO^r)+ra{K^) 



(4.20) 



(4.16) 



where Ta is analytic. Although we presently have no 
theoretical support for this "self-energy" scaling form, 
we appeal to the similarity of the denominator to (j4.13p . 
If one matches (14.161) to (|4.18p . one now assumes that 
T(x,y) « 71 a;^ for small x,y while Ta(\i.,uj) ~ kuj'^ for 
small k,uj. The x'^ term yields T w 1)71 | (5 j"^', with v 
given by (|4.1ip as as required. From the analytic term 
we recover a finite n without any scaling constraint on 
z. Standard static scaling, without any unusual analytic 
corrections, is recovered for to = 0. 

To summarize key the results of this subsection, the 
leading behavior of k is governed by the analytic part 
of the free energy, and is unconstrained (by any argu- 
ment so far presented) by the dynamical exponent z.— 
The critical behavior of k is governed by the exponent 
a, and yields only subleading corrections to the analytic 
behavior. 



where 7 is the susceptibility exponent. At small fc,Ci; in 
the superfluid phase, the dynamics is governed by the 
hydrodynamic Lagrangian (j3.29p . Since this is now a 
bulk Lagrangian, with p the actual bulk density, with 
27r-periodic boundary conditions on the phase 0, the 
ipdr4> term integrates to 2T:inpV = 2-KiN, and there- 
fore drops out of e"^"". Recalling that at hydrodynamic 
scales, if} = ipQC^'^ where ipo is the order parameter, the 
resulting Gaussian Lagrangian yields 



G(x,' 



g-i{[0(x,r)-0(O,O)]^> _ g-([0(O,O)]^> 



|^oP(</'(x,r)^(O,0)), 



(4.17) 



in which fluctuations about the ordered state are assumed 
small. In Fourier space one therefore obtains 



G(k,w) 



1^0 I 



Tfc2 



(4.18) 



B. Crossover exponent associated with 
particle-hole symmetry breaking 

Since we expect the presence of ^(x) to change the 
universality class of the the phase transition, there must 
be an associated positive crossover exponent, </)g, which 
quantifies the instability of the classical random rod fixed 
point with respect to this term. What is the value of (t>g, 
and what conditions does it place on the values of the 
classical fixed point exponents? 

To begin to answer this question, let us write Cc — 
Co 



Cg, where 



d'^X 



dr 



(4.21) 



We assume (see below) that [g(x)]av = 0, [g{x)g{x')]s_v = 
Ag(y3(x — x'), which implies a statistical particle- hole 
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symmetry, and that g{'x.) and r(x) are statistically in- 
dependent. The correlation function takes the form 
ip{x.) = (5(x) for uncorrelated disorder, but for reasons 
that will become evident below, we shall allow more gen- 
eral long-range power-law correlated disorder with 



(^(x) 



'{d+a) 



(4.22) 



for large |x| and some exponent a. The crossover ex- 
ponent, (j)g, which in general will be a function of a, is 
defined by the following scaling form, valid for small Ag : 



fs ~ A\d, 



BAg 



(4.23) 



where ao is the random rod (quantum) specific heat ex- 
ponent, and the subscript on Sg is to serve as a reminder 
that Ag will also generate a shift in the position of the 



critical point: 



Sg=5 



The value of (/)g may 



now be inferred from the derivative 



fdfs 



9/ A„=0 



A|<5|2-"o[S|^|-^.. -(2-ao)ci|<5ri], 
\i\ ^ 0, (4.24) 



where we choose A and B so that <i>g(0) — $g(0) = 1. 
Note the very singular term generated by the shift, 
which may often dominate the 1(51"*^" term of inter- 
est. Now, this derivative may also be calculated directly 
within perturbation theory: 



1 



/(A,)-/(0) = - — 



(^.>o + i[(4>o-(^.>o]+0(5^) 

(4.25) 

where the averages are with respect to £0. Assuming 
that (jf(x) and r(x) self average, we obtain 



/(A,)-/(0) ^ - — 



[(/:s)0]av+2[(4>0]av+O(/) 



(4.26) 



= --A,^(O)[(|V'(O,0)P)o]^^--A, / dr / dV(x)[(V'*5.V'(x,T)V'*9.V(0,0))o]^, + 0(A2), 



where independence of .g(x) and r(x) has been used. 
Thus 



where for 5 > we have lim„ 



-*idy(w;S) — while for 



( df 



9 / A„=0 



= (^(0)£o+ ^ dT ^ d'^xip{x)gg{^, t) 



(4.27) 

)nl , and we have defined the correla- 



\9A, 

where eq 
tion function 

Ggi^, t) = UrdM^, T)rdrHO, O))o]av. (4.28) 

Let us define the Fourier transforms 

GgiKuj) = JdrJ d^xe^'^'^-^+^^^Ggi-K^T). (4.29) 

Then when ^(x) = we have from (j4.6p . 

=eo + ^9(0,0) ^ |(5|(''--°)'^°, (4.30) 

where the exponents are appropriate to the random rod 
problem. More generally we expect a scaling form for 
small |k| and lo: 

T,(k,c.) = £0 + ^5(k,c.) « A^\5\'^''-^«>°y{k^,u^^) 

(4.31) 



(5 < we have \iviiu}^s^oy{'w, s) = 1. Thus, 

= ^^(k)T,(k,0), (4.32) 

where we have used the convenient shorthand notation 
/k = / (Cf' L = / I^' etc. For uncorrelated disor- 
der, </3(k) = 1, while for power law correlated behavior, 
(^(k) w Cafe" -f Co -l- . . . as A: — > 0, where a was defined in 
(|4.22p . The final fc- integral must be treated carefully to 
extract its J-dependence. Let us rewrite 



y{w) 



yiw, 0), yiw ^ CX3) yo > 0, (4.33) 



where the finite limit follows from the requirement that 
T(k, 0) should be a well defined function of k alone (in 
this case a power law oc yok'^~^) when i5 = 0. Clearly 
one cannot simply scale ^ out of the integral in (|4.15p 
since the resulting integral over (p{w/£,)w'^~^y±{'w) will 
not converge. Rather, one must first subtract out the 
large w behavior. Let us write 

yH « yo + yiw-'^' + y2W~^' + yaw-'^^ +... (4.34) 

where the source of the spectrum of exponents, < wi < 
UJ2 < < will be discussed below. Each term 



23 



yjW yields a contribution 

b, EE ^-'^'y, f ^(k)fc-^+^-^o, (4.35) 

SO long as Uj < Wmaxio-) = 2d — zq + min{a,0} (i.e., 
the integral converges at small k; convergence at large k 
being ensured by the implicit lattice cutoff, k < k^ ^ 
Let us define the subtracted scaling function. 



sy{w) = y{w) ~ ^ y-jw' 



Then 



^{\<i)y{kOk 



d-ZQ 



(4.36) 



(4.37) 



+ / fc'^-^«^(k)<5i^(fcC), 
where the last term may be evaluated as 

fc''-"°^(k)(53^(A:0 ~ 6a|<5p''~"°+"^" + 6o|5|^''*""''^", 

ba = Ca [ W^-'o+agy^^^ 
J w 

bo = CO f w^-'-^Syiw). (4.38) 



In fact 6o = in our case because, tracing through 
the definitions, the integral is dominated by the large 
r limit of Gg{0, t), which is proportional to the square of 
{ip*dr^)o- The latter vanishes by particle-hole symmetry 
of the random rod model. 

Now, the origin of the exponents ujj is as follows: the 
operator 



d'^x dr dr'V'*9rV'(x,T)V'*a^V'(x,T') (4.39) 



P 



will have an operator product expansion in terms of 
eigenoperators of a renormalization group transforma- 
tion near the critical fixed point of interest: P = 
hiOi + /12O2 + h^O^ -|- . . ., and O,; is assumed to have 
renormalization group eigenvalue A^. This implies that 
(P) - hi\5\'-'^+''"-^'> + h2\S\'-'^+^°-^^^'' + .... But since 
(P) = J^^y{k^)k'^-'^°, comparison with (1434]) and (|4:35l) 
implies that ujj ^ d + zo — Xj . The loj therefore reflect 
the renormalization group transformation properties of 
P near the fixed point. 

We can now understand the crossover exponent, (j)g. 
Comparing ([iTM)) and ^j^^ together with (with 
bo — 0), we see that 



2-ao 



V min{a>j ,2d — zo + a} 

max{Aj, 2^0 - d - a}. (4.40) 



The crossover exponent is therefore either 4>g = Xiv, or 
for small enough a, (f>g = {2zo — d — a)v. Thus 



implying 



2zo 
Ai, 



a < 2zo — d 
a > 2zo — d - 



Ai 
Ai, 



> ^ 



zo > 
Ai >0, 



a < 2zo — d — Ai 
a > 2zo — d — Xi. 



(4.41) 



(4.42) 



In particular, for short range correlated disorder, where 
in effect a ^ 00, we require Ai > in order that dirty 
boson disorder destabilize the random rod fixed point. 
The exponent Ai is a nontrivial exponent, and we shall 
compute it within the e, er-expansionJ^ 

A naive estimate for Ai is obtained by supposing that 
the equality 2zo — d — a = Xi should occur when a ~ 0, 
yielding Ai ~ 22:0 — d which becomes positive for zq > f ■ 
Note that this same estimate would have been obtained 
from the second term in (|4.2ip if we had assumed bo ^ 0. 
As an aside, this estimate is actually exact in the cor- 
responding derivation of the Harris criterion for classi- 
cal disordered magnets. There the correlation function 
(|'!/'(x)p|i/)(y)P)o appears. Since (ip^)o does not vanish, 
neither does the coefficient analogous to bo- This gives 
rise to a free energy contribution (|'0P)o ~ J^"^"" which 
leads immediately to the Harris criterion, (j)g = ao- 

In fact, we shall find that Ai > 2zo — d, i.e., a drops 
out at some negative value, and (j)g becomes positive for 
Zo larger than a value somewhat less than |. The ran- 
dom rod result zq = 1 in = 1 is consistent with this 
criterion, although this case is somewhat special because 
the random rod fixed point is the same as the pure fixed 
point in d = 1 (i.e., random rod disorder is irrelevant, 
though boson disorder is relevant). The generalized Har- 
ris criterioni^ indicates that rod disorder is irrelevant 
when a 



pure 



^purc 



< 0. Using hyperscaling (valid here 
for d < 3), and the fact that z = 1 at the pure fixed 
point, this requires impure > | (compare the less strin- 
gent requirement, t'pure > with dtot = d -I- 1, for the 
usual Harris criterion for point disorder). For d = 1, 
t'puro — > 00, while for d = 2, impure — |, so the pure fixed 
point becomes unstable to rod disorder somewhere in be- 
tween d = 1 and d — 2. In all cases where rod disorder is 
irrelevant, one then trivially has zq > |, and dirty boson 
disorder will certainly be relevant. 

We now turn to the question of the relevance of go = 
[g(x)]av, i.e., of full breaking of particle-hole symme- 
try. If one carries through a naive scaling analysis using 
^90 =Id''xf dr (iffglVP - gordri') in place of gUJ, 
one obtains df/dgo = 0, while 



-^=eo + gg{0,0) = Tr^\5\ 
(^9o 



(rf-^o)i'o 



(4.43) 



If the go-dependence of the singular part of the free en- 
ergy scales in the form 



50 



(4.44) 
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it is tempting to identify the crossover exponent via 

2 - ao - 2(f)gg = {d- zo)v (l)go = zqv (4.45) 

which is always strongly positive. In the absence of all 
disorder, this is the correct exponent (with Zpurc — 1, and 
^'puro = '^d+i) describing the crossover from the {d + 1)- 
dimensional XY behavior at the tip of the Mott lobe, to 
the generic onset of superfluidity in a dilute Bose gas as 
the density is increased from zero away from the tip. This 
same exponent describes the shape of the pure system 
Mott lobes near their tips. Fig. [Ija): since the transition 
line is defined by some critical value Xc of the scaling 
function argument x — go/|i5|'^"', this leads to go c ~ 
Xc|<5|^«o, i.e.^ 



ficiJo) - \Jo - J°\ 



XY 
'd+1 



The same argument implies that 



(4.46) 



(4.47) 



near the random rod critical points in Fig. [TJc). It is 
likely that zqVq > I in d = 3, leading to the pictured 
cusps, 

One must be careful in interpreting the result ()4.45|) 
in the presence of random rod disorder. The random 
rod problem leads to an incompressible glassy phase [see 
Fig. [He)] . Perturbing the random rod critical point with 
either or Ag therefore not only changes the critical 
behavior, but also the nature of the glassy phase. The 
difference between (j)g in (|4.40p and 4>gg in (|4.45p therefore 
reflects the different rates at which the crossover to the 
Bose glass phase occurs under the influence of the two 
perturbations, with 170 clearly having the stronger effect. 

We emphasize that the two crossover exponents de- 
scribe, in renormalization group language, the rate at 
which one is initially driven away from the random rod 
fixed point in the two orthogonal directions described by 
go and Ag. They tell one nothing about the eventual 
termination of the flows on some new stable flxed point. 
We have argued physically, supported by the droplet pic- 
ture of Sec. mil that although initially growing rapidly, 
go must vanish again as the dirty boson flxed point is 
approached, its main role being to induce a finite value 
of Ag under renormalization. Quantifying this expected 
irrelevance of 170 (i.e., 4>g„ < 0) at the commensurate 
dirty boson fixed point )^ would require an analysis of 
the right hand side of (|4.43p in the presence of a finite 
value of Ag. We saw in the previous subsection that in 
this case = k is dominated by analytic terms in the 
free energy, and hence that (jjgg is related to subleading 
terms in the singular part of the free energy that cannot 
be inferred from a simple scaling analysis. In Sec. I VII will 
obtain ^g^ within the e, er-expansion, and demonstrate 
explicitly the manner in which particle-hole asymmetry 
becomes irrelevant at the commensurate fixed point for 
sufficiently large Cr. 



V. CALCULATIONS IN ONE DIMENSION 

A. Sine-Gordon model 

In this section we review and expand upon the analy- 
sis of one-dimensional versions of the dirty boson prob- 
lem, with emphasis on the weak disorder limit i^i^ In 
App. [B] we derive various dual representations for the 
one-dimensional Lagrangian based on the discrete-time 
Villain representation, (|B4p . We shall analyze the sine- 
Gordon version, (|B14p with (jBlSP : 
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^M/(9/^R) + 2yo^cos(2^5R), (5.1) 
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where we have assumed that Vjj = VqSij is diagonal. 
Here R = {I,T), with integer I,T, are points on a dis- 
crete space-time (dual) lattice, — 00 < S'r < 00 are con- 
tinuous spin variables, the discrete derivatives are defined 

by 9/ S'r = S'(/+i_T) — 'S'(/,T), i9t<S'r = S^j^t+i) — 'S'(/,t), 
and the cosine term represents an external periodic po- 
tential which prefers integer values of S'r. This model 
has the physical interpretation of a fluctuating interface, 
represented by the "height" variables 5r. The coeffi- 
cient Kj is proportional to the Josephson coupling [see 
equation (|Bip ]. In the absence of /z/, which has the inter- 
pretation of a random tilt potential, the phase transition 
in this model is from a flat phase at large 1/ Kj, where Sr 
has only small fluctuations about some integer value and 
exponentially decaying correlations, to a rough phase, at 
large if/, in which the interface wanders and has logarith- 
mically divergent height-height correlations. This rough 
phase corresponds to the superfluid phase in the boson 
model, and the renormalized, long wavelength value of yo 
vanishes. In the presence of the random tilting potential, 
/i/, the rough phase is qualitatively unchanged, but the 
flat phase is no longer necessarily quite so flat: see below. 

Let us decompose /j,/ = ^o + ^l^i into a uniform part /io 
and a random part with [<5/i/]av = 0. Thus, nonzero /io 
represents the breaking of particle-hole symmetry. When 
/io = the interface will be globally flat. For sufficiently 
large /io > /io,c, where /^o.c represents the Mott gap 
(which will vanish for some combination of sufficiently 
large J/x/, large Ki, small Vq and small yo) the interface 
will acquire a global tilt, with [(Sr — S'R')]av — Qo{I~I')- 
The exact form of the slope function 



Qo(MO,^^0,K),yo) = [(5/S'R)]a 
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9/io 



(5.2) 



which is proportional to the density difference from the 
Mott phase, can be computed perturbatively in powers 
of yo (see below). In the special case m/Vo = \, the 
spatial derivative terms in (|5.ip may be combined in the 
form Vb(9/SR — 5)^, and there is an exact degeneracy 
between 9/Sr = 0, 1 (the integer values preferred by the 
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cosine term). This is the particle- hole symmetric model 
at half filling, Qo = ^^'^^ 



B. Perturbation theory in the superfluid phase 



Let us define the random walk 



Wl 



Vo ^ 

' J=0 



(5.3) 



(defined to be minus the sum from J = / to for J < 0) , 
and let 



Sr^ Sr- QqI -wi, R 



With this choice, one will have [(S'r — S'r 
sine-Gordon Lagrangian takes the form 



(/,T). (5.4) 
av = 0. The 
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{drS-ELf + VoidiSnf 
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2j/o^cos[2^(^R 
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wi + Qol)] 



+ Eo(3L, 
where the constant term is 



(5.5) 



{QoVo + S^iI){2flo - QoVo + Sfii) 



R 



21^0 



2Vo 



(5.6) 



In (j5.4p we have dropped a sub-extensive boundary term 
-(/xo - QoVo)J2R9iSn since Yl^ndiSn = J^ri^L.T - 
So,t) — 0{\/l3L), and therefore yields vanishing contri- 
bution in the thermodynamic limit. This term may in 
fact be made to vanish identically by choosing periodic 
boundary conditions for S'r. 

The explicit dependence on /io is only in the last (con- 
stant) term of (|5.5p . and if follows from (|5.2p that Qq 
may be determined determined by minimizing the free 
energy at fixed /io: 
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(5.7) 



When 2/0 = the condition (|5.7p yields Qo = A*o/K) and 
the SjjLi yield only a trivial additive constant to the free 
energy. In this limit, for Kj = Kq fixed, the two-point 
correlation function is given by 



G(R-R') = -((5r-SrO 



(5.. 



In 



p(R - R') 

Po 



oo, 



where, 



p(R-R') = 



KoVo 



(5.9) 



is the appropriately rescaled distance, and po = 0(1) is 
a constant scale factor. When Kj fluctuates, its disorder 
average must be included. The result is still (|5.8p . but 
Kq then becomes a complicated effective parameter. The 
generalized Harris criterion (Ref. 15 and Sec. lIVp implies 
that disorder in the coefficient Kq is an irrelevant per- 
turbation at the pure (2D XY) critical point at integer 
filling in d = 1, so we will, for the rest of this section, 
simply take Kj = Kq when considering the influence of 
the 

Let us then consider the uq term as a perturbation on 
the quadratic term in Csg- Deep in the superfluid/rough 
phase, where Kq/Vq is large, this is a well defined expan- 
sion. It is also well defined when Qo is not too small: 
the cosine term in (j5.4p then oscillates very rapidly from 
site to site, and effectively averages itself out. This cor- 
responds to the region between Mott lobes in Fig.[T] The 
condition (15.711 leads to 



Qo = j^- CQ{fio,KQ, Vq, A)yl + 0(2/4), (5.10) 

where the positive coefficient of the correction term is 
given by 

A-TT 

Co^ttY.^ [(sin[27r(SR ~Sq + wi + QqI)])q],, , 
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(5.11) 

in which the thermodynamic average is with respect to 
the Gaussian Lagrangian £sg(2/o = 0). To evaluate this 
further we assume that the 5pi are independent, with a 
symmetric distribution. Let us define a measure of the 
disorder strength. A, via 



e ^0 ^' 



(5.12) 



Then, 



Co ^ Ve-4-^«We-2'^^^^l^l/sin(2^/.o//^o). 

(5.13) 
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The sum over I clearly converges. Using (|5.8p . it eas- 
ily s een tha t the sum over T converges so long as ujq = 
2TTy/ Kq/Vq > 1. We shall see below that the superfluid 
phase is define by wo > 3, so this condition is indeed met. 

The corrugation due to t/o therefore slows the rate of 
climb of the interface from its unperturbed rate, p/Vq. 
When Kq/Vq and fiQ become small this perturbation the- 
ory breaks down — a signal of the phase transition into the 
Bose glass phase. 



C. Stability of the superfluid phase 

We consider now the stability of the superfluid phase to 
Uq. In Ref. this analysis was performed using Kosterlitz- 
Thouless-type renormalization group methods. Here we 
will take a less sophisticated route and adapt the scaling 
approach described in Sec. IIVBI This calculation will 
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also allow us to examine the effects of finite /iq on this 
stability. The latter will allow an explicit confirmation of 
the irrelevance of full particle-hole symmetry breaking at 
the commensurate (statistically particle-hole symmetric) 
critical point. 

We consider the relevance of the cosine term on the 
fixed line, characterized by the long-range correlations 
(15. 8p . To this end, define the local operator 



Or = COS 



27r(^R + WI + Qol) 



(5.14) 



and introduce a "temperature" variable, analogous to S 
in Sec. llVi by adding a mass term 



(5.15) 
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to CsG- By this device we may discuss the relevance of 
the yo term to the critical behavior as t — > 0. To this 
end, we postulate a scaling form for the singular part of 
the free energy. 



fsiyo) ~ At 
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so that 



1 (d^fs\ 



'$'(0). 



(5.16) 



(5.17) 



The superfluid phase always occurs at = 0, so there 
will be no shift in the critical value t — 0. As usual, the 
2/0 term is relevant if (j)y > 0. 

The derivative in (j5.17p may be computed in terms of 
the average 
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cos(27r(5o-f), 
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where (I5.12p has been used, and where 



G(R,t) = 



n -1 



R,0 
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For i ^ 0, G(R, t) has the logarithmic form ([5 
finite t one may write 

g-4.=G(R,t) ^ [p(R)/p„]-o^[^(R)2^/^2j^ 



(5.19) 
For 



(5.20) 



where = 27r-^ Kq/Vq determines the power law decay 
of correlations at criticality (i.e., in the superfluid phase) 
and the scaling function E{w) decays exponentially for 
large w [this can be seen explicitly by writing G(R, t) = 



G(R, 0) + 5G(R, t) and using ((5T9)) ] and E{Q) = 1. This 
exhibits the scaling of the correlations when ?/o = Oj Etnd 
and since p scales with \/t one immediately identifies the 
correlation length exponent ly = ^. 

In addition to the subleading singular part we seek, 
(|5.18p contains analytic terms in t, whose Taylor coeffi- 
cients may be evaluated by taking derivatives of ()5.18p 
with respect to t at t = 0. Let n be the first positive 
integer such that the nth derivative of (|5.18p diverges 
as t 0. Since each derivative of (|5.20p with respect 
to t brings a factor of /5(R)^ ^ |Rp out of the scaling 
function, this divergence arises from a failure of the inte- 
gral to converge at infinity. The leading singularity may 
therefore be computed exactly by considering only the 
large |R| asymptotic behavior of the integrand. In this 
limit one may perform the strongly convergent sum over 
/ by setting 7 = inside G. Defining, 



Do = -8^cos(2^go/)e-'"'^'l'l, 



(5.21) 
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one finds 

X E^^HKoVoTH/pD 
« EoDo&"-^^l^~'^, (5.22) 

where E^^\x) is the nth derivative of E{x), Ti — 0(1) is 
a lower cutoff [whose arbitrariness yields only subleading 
corrections to the last line of (|5.22p ] , and with coefficient 
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It is clear at this point that the n > we seek is the first 
integer for which cjo ^ 2n — 1 < 0. One finally obtains 
the singular part 



dVo 



,o=0.sing r[(^o + l)/2] 

(5.24) 

From (|5.7p and (|5.8p we see that, up to scale factors, 
space-time is isotropic. Thus z = 1 and hyperscaling 
yields 2 — a = 2iy, so that from (|5.17p we may finally 
identify 

0y = = 1- ^VKo/Vo. (5.25) 



Hence yo becomes relevant when y/Ko/Vo < This 
should be compared to the analogous result, y/Ko/Vo < 
-, for the usual Kosterlitz-Thouless transition where 
Pj = 0. Thus, the interface roughens earlier (i.e., at 
smaller Ko), meaning that superfiuidity is more stable, 
in the presence of disorder. For y/Ko/Vo > is 
irrelevant and may be set to zero to calculate universal 
quantities near the phase transition. At the critical point 



27 



one has uj = lUc = 5, which should be compared to the 
Kosterlitz-Thouless value, uic = 4. One may then, for 
example, invert the duality transformation in this limit 
to obtain the actual superfluid correlation function. One 
finds that jU, with (|B16p . takes the form 



Vo 



0{<Pr+i 



Ur+f 



(5.26) 



where r = (i,T) is the direct lattice integer position vec- 
tor and where now — oo < (/)r < oo is a continuous phase 
variable [since (jB16p forces Vxm=Oas?/^0, we may 
write m = Vp, where p is an integer scalar field, then 
define ipr — fpr ^ '^T^Pr] ■ Thus 

(5.27) 

where p(r) is the same as p(r) in (|5.6p . but with XoVb re- 
placed by j^^y^ ■ The exponent rj is defined in such a way 
that g{i,T = 0) - \i\-{d+z-2+-n) criticality. Equation 
(|5.27p then follows since d — z = \ and p(i,T = 0) cx |i| 
for large At the critical point we have 77 = i, which 
should be compared to the usual Kosterlitz-Thouless 
value, V = J- 

The above calculation was performed at j/o = 0. When 
yo > 0, in the region where it is irrelevant, the parameters 
Kq and Vq in the Lagrangians (|5.5p and (|5.26p must be 
renormalized to values -R'aXj/o) ^^nd Vr(?/o) before setting 
?/o = in the derivation of (|5.27p . Thus oj = 2tv ^JKr/Vr 
and 0j/ = I — tta/ Kh/Vr, but the relation 77 = i is 
still exact. The parameters Kj^ and Vr are the exact, 
long wavelength (hydrodynamic) interface stiffness mod- 
uli that a bulk experimental probe would measure, and 
are directly analogous to the superfluid density and com- 
pressibility in the superfluid problem — see (|3.29[) . The 
above analysis shows that when the ratio \JVrIK-h ex- 
ceeds the universal value yo becomes relevant, and 
simple renormalization of the Gaussian Lagrangian (|5.26p 
is invalid. We then expect Vr/ Kr — > cxd, and the inter- 
face becomes localized. At the critical point separating 
the localized and delocalized phases, the interface is still 
delocalized, with the universal parameter values quoted 
above. Using renormalization group techniques, all of 
these results may be confirmed by constructing the de- 
tailed flows around this fixed point 1^1^ 



D. Restoration of particle- hole symmetry in ID 



Recall now the discussion in Sees. II CI and IIVBI 
of asymptotic restoration of statistical particle-hole 
symmetry — namely the irrelevance of /xq in the presence 
of nonzero This is seen trivially in the ID case be- 
cause the critical fixed point occurs at y^ = 0, at which 
point the mapping 5r = S'r, Eq. (|5.4p with Qq — fiQ, 



entirely eliminates ^0, as well as all the 5^/, from (|5.5p . 
except for the analytic additive term (|5.6p . One there- 
fore obtains in this case a rather extreme form of irrele- 
vance, in which the influence of /io does not decay with a 
characteristic exponent (l)g„ < 0, but actually disappears 
entirely. 

More generally, when both /ip and yo are nonzero, the 
fact that /iQ appears only in the cosine term in (j5.5p 
means that its influence must vanish on large length 
scales whenever j/o is irrelevant. Examining the scaling 
analysis (|5.17p - (|5.25p . used to determine the range of this 
irrelevance, one sees that ^0 (via Qo) appears only in the 
cosine factor in (|5.2ip . This factor is completely dom- 
inated by the exponential decay due to the fluctuating 
part of the ///, and is therefore of no real consequence, 
producing only analytic corrections multiplying the lead- 
ing singularity, and therefore having no influence on the 
value of the crossover exponent (j)y. The ultimate ori- 
gin of this result can be seen in (|5.14p : only the value 
of 9 1 = wj mod 2tt is important, and when the vari- 
ance measure. A, of (5/i/ is sufficiently large this field is 
basically uniformly distributed over the interval [0,27r), 
irrespective of the "mean drift" iiI/Vq. 

Note, finally, that the irrelevance of a uniform /zq has 
no bearing on the value the dynamical exponent z, and 
it can easily be verified that the finite piece in the com- 
pressibility comes purely from the analytic part of the 
free energy, especially the term (|5.6p . Thus, although 
z — d — 1 in this case, one may view this as a 'coin- 
cidence' that has no bearing on the general mechanism, 
discussed in detail in Sec. IIV A 21 by which k remains 
finite through the Bose glass-superfluid transition. 



VI. THE EPSILON EXPANSION 

In this final section we turn from exact calculations 
in one dimension to approximate calculations in higher 
dimensions, expanding on, and providing more context 
for, our previous worki^ Unlike the classical point disor- 
der problem, the classical random rod problem [ (j2.4p . or 
(|4.ip with g{x.) = but 7-(x) random], does not have a 
simple epsilon expansion about d = 4. Rather, as shown 
in Refs. Ua, one must consider also the limit in which the 
dimension, er of the rods is small, and perform a double 
expansion in e = 4 — D and er (recall that D = d + e^- 
is the total dimensionality). The exponents take mean- 
field values, z = 1, = i, 7/ = 0, etc., at e = = 0, 
and deviations from these values may be computed as 
two- variable power series in e and er- 

Our purpose in this section is to extend this technique 
to the dirty boson problem. We saw in Sec. IIVBI that 
a certain nontrivial crossover exponent 4>g must be pos- 
itive if, as expected, particle-hole symmetric disorder is 
to lead to new critical behavior, different from that of 
the classical random rod problem. This result was con- 
firmed explicitly for d = 1 in Sec. [V] there, random rod 
disorder was found to be an irrelevant perturbation on 
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the pure (Kosterlitz-Thouless, 2D XY) critical behavior, 
whereas dirty boson-type disorder was found to be rel- 
evant, leading to new critical behavior. We shall find 
that for small particle- hole symmetric disorder is an 
irrelevant perturbation on the random rod problem, and 
therefore that the crossover exponent changes sign, from 
negative to positive, at a certain value, e-r = e^(_D). To 
first order in er we obtain the estimate e^(Z) = 4) = ^ 
{D = A yielding d = 3 at e,- = 1). For Cr > e-% there 
are then two fixed points, the stable dirty boson fixed, 
and the unstable random rod fixed point. This then es- 
tablishes the nonperturbative nature of the dirty boson 
fixed point. 



A. Scaling for general e-r 

Let us begin by extending the scaling arguments to 
noninteger e,-. We consider the following generalization 
of (l4T|) (or, equivalently, of £5 in Table ij): 



(6.1) 



where g(x) is an eT-dimensional vector. This form is 
based on (|2.9p . with the same simplifications used in 
(|4.1[) . We write g(x) = go -I- i5g(x), and assume that 
(5g(x) is isotropically distributed in r space. This yields 
the correct er = 1 limit, and ensures that the free en- 
ergy depends only on go = |go|- Clearly, go = is the 
generalization of statistically particle-hole symmetric dis- 
order. As before, we also write r(x) = ro + 5r(x), with 
[Sr]iiv — 0, and ro is the control parameter. 



1. Stiffness constants and hydrodynamic action 



The evaluation of the stiffness constants, (|3.14|) and 
(|4.6p . is slightly more complicated now: although the spa- 
tial stiffnesses, T^, are as before, the temporal stiffness 
now takes on a tensor character. Consider a 0-boundary 
condition in the x-subspace: 



^(x, T + f3Tf,) = e'*'^ V(x, -r), = I, . 
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(6.2) 



Defining the periodic field, ijj = e ^^''^^^ip, and substi- 
tuting into (|6.ip . we find that 

£,[V';go] = /:c[^;go- W]- (6.3) 

The free energy, f = -/J-^-L^'^ In {tr }, is 

therefore shifted by 
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(6.4) 
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Isotropy implies that 
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where we have defined the scalar quantities 



Po 



and (|6.4p reduces to 
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It is straightforward to write down expressions for po 
and Kq analogous to (|4.6p . but we will require only 



Ko(go = 0)^T,= [d^nj^^ (6.8) 

d^X I d'^T[{rdrM^,T)rdrMO,0))], 



where ri is any given direction in T-space. The long 
wavelength action generalizing (|3.29p then takes the form 



•S'cff — — 



|^[|V.0r-(go-V.</.)^] 
2.90 

^Ko(go ■ V^0)2 -I- ipo(go ■ ^t(, 



(6.9) 



Note that when go we have po/ffo ~^ '*o(3o = 0) = 
T,-, and ^off reduces to the more familiar form 

5'cff(go =0) = -^Jd''xJ d'^T (T|V0p + T,|V,0p) . 

(6.10) 

Defining the frequency variables t^n = gQ ■ u: and ujj_ = 
UJ — W||gOi and following the derivation of (|4.18p . equa- 
tion (|6.9p yields a long wavelength, low frequency Green 
function 



G(k,c.) ^ (|^(k,a;)|^ 
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(po/.9o)|'^±P + «oc^^ + T|k|2 



, (6.11) 



where ^po = [{ip{x,T))]^^ ~ |^|'^ is the order parameter, 
and for slow variations, ipi^i ''") ~ V'oe*'^^''''^-'. In deriving 
this form we have neglected the surface term, ipogo • V^-^, 
relative to the others. This is valid in the superfluid phase 
where (j) has only small fluctuations about the long range 
ordered value, which we have taken to be = 0. As we 
saw in Sec. lIII C^ in the disordered Bose glass phase this 
term does become important (see also below). 
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2. Dynamical exponents and scab 



Equation (|6.1ip determines the correlations in the hy- 
drodynamic Umit, even near the critical point. One could 
propose a scaling function of the form: 



G(k,a;) « Goe'-'',g(fce,c.ier ,c.||e^), 



(6.12) 



where k, ujj_ and appear scaled by their appropriate 

correlation lengths, ^, ^| ~ ^^ii and ~ f^-^, where, 
for go 7^ 0, we have allowed for different scalings parallel 
and perpendicular to go- When go = we expect ^| = 
= and z\ \ = z± = z. To be consistent with (|6.1ip . 
equation (I6.12|) requires that 



-, x,y,z^ 0, (6.13) 



in which g^, gi_ and are universal numbers. Thus 
no « Go-VoP.9||(d)'r'+" 

po/go « Go-VoiV(e.^)'r'+'' 



(6.14) 



[compare (|4.18p and (|4.19] . The hyperscaling relation is 
now 2 — a = [d + Z|| + (e^ — l)zj_]i/. Along with the 
usual scaling relations, a + 2/3 + 7 = 2 and 7 = (2 — 77)1^, 
this immediately implies that T ~ \S\^ , kq ^ \S\^^^ and 

P0/.90 ^ l^r^ : where S = ro — ro,c(5o) is the deviation 
from criticality, and the implied exponent relations 



[d+{er - l)z± + Z|| - 2]v 

[d+ {er - l)z± - 

[d + (er - 'i)z±_ + z\\\v, 



(6.15) 



generalize (|4.11l) . If z±_ ~ z\\ ~ z, then v — {d + ~ 2)v 

and vi = = Vt — [d + {er — 2)z\v. If kq remains finite 
through the transition even for noninteger e,-, then we 
would predict 



(?) 



(6.16) 



This is the natural generalization of the proposed scaling 
relation z = d aX €t = liSiSS 

However, as discussed in Ref. |2^ and Sec. IIV A21 this 
argument actually fails for e,- = 1, and there is no reason 
to expect it to fare any better for general e,-- In partic- 
ular, the absorption of the twisted boundary condition 
into the simple shift, (|6.3p and (|6.4p . of go again implies 
that, at least for g^ ^ 0, the critical behavior of po and 
kq is not tied, as in (|6.12p . to the correlation times, but 
only to derivatives with respect to g^. Thus, although 
the relation for v in (|6.15p is expected to be correct, the 

relations for vi and are not. Rather, as discussed 
above (|4.14p , the leading singular contribution to po and 



kq will come from the go dependence of ro^dgo), leading 
to 



(6.17) 



We will see that the renormalization group analysis gen- 
erates an independent value for z, disagreeing, in partic- 
ular, with (|6.16p . confirming again that there is no \S\^^ 
contribution to the compressibility. 



3. Relevance of particle-hole symmetry breaking 



We next calculate the crossover exponent (pg associated 
with nonzero ^g at go = 0, for general e-r- The calcula- 
tion is essentially identical to that leading to (|4.40p . The 
perturbation, (|4.2ip . now becomes 

C^ = J d^x J d^-r|i|g(x)nV;|2-V'*[g(x).V.]v|. 

(6.18) 

We assume [g^{x)g^{x)]^^ = Ag(?!)(x - x')(5pi., where (j) 
is a delta function for short range correlated disorder, 
and varies as for large x for power law correlated 

disorder. Equation (|4.32p now becomes 



f df 



(^(k)T,(k,0) 



(6.19) 



in which T^(k,a;) is defined by (|i:^ - (|i3T|) . but, as 
in (|6.8p . with time derivatives dr in (|4.28p . and 

random rod compressibility exponent ?Jr,o(eT) — [d+{eT — 
2)z]h'o replacing (d — zoji^o in the scaling form (|4.3ip . 

As in (|4.33p and (|4.34p . y{w) will have a similar spec- 
trum of exponents, < a;i(e^) < 072(61-) < i^si^r) < ■ ■ , 
now depending on e,-. This leads in the same way 
to and (133H1), still with 60 ee 0, and a term 

^^|j|[2d-i-(e,-2)2o]f. We then find 



(I>g/VQ 



2zq ~ d~ a, a < 2zq — d — Ai(eT-) 



Ai(eT), 



a > 2zo — d — Xiicr) 



(6.20) 



essentially as before [see (|4.41l) ]. but now with Xj = d + 
StZ — Uj and all exponents evaluated at the random rod 
fixed point in er time dimensions. 

We shall compute Xi{er) to 0{er) below. From the 
naive (and, as we shall see, incorrect) estimate, Ai = 
2zo — d, we again have Ai > for z > |. Since z = 1 
at er = 0, we expect, as stated earlier, Ai < for small 
er, becoming positive only for e,- > > 0. We shall find 
that for er > e'^ a, new stable fixed point with A* > 
bifurcates away from the random rod fixed point (with 
A* = 0). The exponent zo is substantially smaller than 
I at this point, violating (|6.16p for any er > 0. Assuming 
that this new fixed point may indeed be identified with 
the true dirty boson fixed point when = 1, we conclude 
again that l|6.16p is incorrect. Below we shall generalize 
the analysis of the excitation spectrum of the Bose glass 
phase in Sec. |lIl]to general er- For er < 1, the issue of 
whether the Bose glass phase is compressible turns out 
to be rather subtle. 
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B. Phase diagrams for general er for any set of integers n = {UaY^^ 



To understand the phase diagram for general e,- it is 
instructive to first consider the mean field approximation, 
in the form of the infinite range hopping model {Jij = 
Jo/N for all where N is the number of sites). In the 
absence of disorder, and with onsite repulsion only, the 
mean-field approximation to the Lagrangian (j2.7p takes 
the form 

Cmf = J d'^rl ^^cos[0,(T)-0j (r)] 

- ^E(^-'^' -*/^)'}' (6.21) 

with temporal periodic boundary conditions 

0,(T + /3e„) = 0,(r) + 27rnc„ (6.22) 



1. Kac-Hubbard-Stratanovich transformation 

A complex Kac-Hubbard-Stratanovich (KHS) field 
M{t) may now be introduced to decouple the hopping 
term. Using the identity 

^a\z\-/N = d2^g-a(^.*+V.*.+Ar|^|^)^ (6.23) 

where the integral is over the real and imaginary parts of 
■0, the partition function can be written in the form 

ZMF = y"^M(-r)e^^"^[^l, (6.24) 

in which the integral is over all Cr-dimensional complex 
functions M{t), weighted by the effective action 



Jo 
2 



(6.25) 



in which the decoupling has reduced the phase variable 
trace to one over a single variable (/)(t), which includes 
a sum over all possible boundary conditions (|6.22p . The 
constant a = ^Jdf'^r incorporates the discretization of 
the T-integral, and the measures DM{t), D(I){t) incor- 
porate any normalization factors needed to make sense 
of the continuum limit. 



2. Saddle point evaluation 

In the thermodynamic limit, A'^ — > oo, a saddle point 
evaluation of (|6.24p becomes exact. Since the superfluid 
order parameter is homogeneous and static, the lowest- 
energy saddle point must be a time independent field, 
M(t) = Mq, which can be chosen real. Near the criti- 
cal point Mq will be small, and 5mf[Mo] has a Landau 
expansion, 

/MF(Afo) = -r'-SuFlMo] 

= /o(/x, Jo) + ^ro(M, Jo)|Afo|' 

+ ^uoifi, Jo)\Mof + . . . . (6.26) 

The connection between Mq and and the order param- 
eter ipo = {jfJ2j^^'^') is obtained by adding a term 
-i / d'-TY,j {h*e"^^ + he-'-t'') to ^mf- Since this term 
is local, the KHS transformation simply yields the re- 
placement Af M = M + h/Jo inside the logarith- 



mic term in (|6.25[) . In terms of M, the first term 
in (p:^ becomes ^Jq\M\'^ = ^Jo\M - /ip, and the 
free energy is obtained by minimizing /MF(Afo, /i) — 
fMF{Mo) ~~ \{h*MQ + hMl) + |/i|V2Jo. One therefore 
obtains Vo — — 2(9/mf/9/i)/i=o — Mq. Therefore AIq is 
in fact the order parameter in this theory. 

Assuming only that uq > 0, the transition line occurs 
at ro = 0, and the superfluid phase corresponds to ro < 
0. To obtain the phase diagram, we therefore need only 
compute To, which is given by 

ro = Jo- Jo J d'-T{cos[(j){T)]cos[(l){0)])g 

= Jo-ijo'Re J d^^T(e^[*M-^(°)l)^ (6.27) 

in which the average is with respect to the Afo- 
independent ^^j^ [Vr</'(T) — i/x]^ term. To further evalu- 
ate ro, we first account for the boundary conditions by 
changing to the periodic variable 

Ht)^4>{t)~^^. (6.28) 

It follows that (/)(x) - 0(0) = 0(x) - 0(0) -I- 27rn • 
Vr0(T) = Vr0(-r) + 27rn//3, and /d^-rV^^ = 0, and 
one therefore obtains 

y- -0'-{2^n/|3-i^,f/2Uo„2■K^n■T/|3 
/ i[0(-r)-0(O)]\ ^ Z^nS f 

\ /o ^^g-/3'-(27rn//3-j,*)V2C/o 

X e-5<W-)-'^(o)]'>oo, (6.29) 
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in which (•)oo is an average with respect to the the 
Gaussian term (Vr'^)^/2C/o, and we have used property 
{e^-^) — e^^"^ vahd for any Gaussian random vari- 
able X. The remaining average is the inverse Fourier 
transform of L/q/I'^P, namely the Green function for the 
Laplacian in Et- dimensions: 



Mean field phase diagram: e = 0.2, 0.6, 0.9, 1 



A{er) 



r(e./2) 



2(2-er)7r'^-/2' 



(6.30) 



Note that ^(1) = \. 

For the n-sums we make use of the identity (13.371) . The 
ratio of n-sums in (|6.29p becomes 
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(6.31) 

in which the right hand side is valid for [3 — > oo, and lo(At) 
is the value of 1 that minimizes (/x -f [/o/3^~'^^lo)^- For 
Er < 1 one clearly obtains Iq = 0, while in the special case 
e-r ^ I one obtains the periodic form Zo(/^) = + 5]- 
Adopting the convention — ^ < 1.1 /Uq mod 1 < ^, we 
finally obtain, 



g-^Uo\T\-UoitJ./Uo mod 1)t — 2. 



(6.32) 

which should be compared to the droplet model calcula- 
tion JSUl), (jOSll . 



3. Superfluid transition line 



Substituting (|6.32|) into (|6.27p . one obtains the follow- 
ing equation for the transition line: 



•^o.c(ImI) 



Jo 



I^i\t^\r), (6.33) 



where Ii^{z) is the modified Bessel function. Note that 
(z/2)~''/y (z) is analytic at the origin, taking the value 
l/r(i^ + 1) there, and that I_i{z) — y/2/TTz cosh(z) re- 
duces (16.331) to the correct result at Cr = 1 • Note that for 
Cr > I the integral (|6.33p diverges for 7^ 0, indicating a 
divergent compressibility. Thus, only e < 1 displays the 
physics of interest to us. 

Using the Taylor expansion of the Bessel function,''"' 



FIG. 8: (Color online) Mean field phase diagram for the 
Josephson junction model, showing Mott (M) and superfluid 
(SF) phases, with particle-hole symmetric transition point (C) 
for e-r — 0.2,0.6,0.9, 1. Larger er correspond to smaller Mott 
lobes. For < 1 (dashed lines) there is a single Mott lobe 
persisting for arbitrarily large fj,. For er — 1 the Mott lobe 
shown (solid line) is actually repeated periodically. The scaled 
axis variables are Jo^dlfJ-l)/ Jo,c{0) and fi, as defined in (|6.34|l . 



one obtains, 

Jo,c(0) 



-2k 



2[UoAier)]- 



[UoA{er)] — 



(6.34) 



It follows that Jo,c(0) — > 2 as Et 0, and Jo,c(0) Uo/2 
as e-r ^ 1. 

For Cr — I the result (|6.33p is valid only for the central 
Mott lobe. The remainder are constructed by periodic 
repetition. For this case, the integral can be computed 
analytically for arbitrary /i, and one obtains;^ 



Jo,c(pt) 
Uo 



1-f^ 



(6.35) 



For Er < 1 the |tP^'^^ term dominates the exponent in 
(|6.33p for large r, and the integral converges for arbitrar- 
ily large A saddle point estimate yields 



exp 



2p. 

2-€r 



Ai>l, (6.36) 



where non-exponential prefactors have been dropped. In 
Fig. [5] we plot numerical solutions of (|6.34p for the phase 
boundary for different values of e,-. 



32 



Since in mean field theory, every site gets effectively 
decoupled, it is very easy to incorporate the effect of site 
disorder within this formalism. The generalization of the 
mean field action is 



is the /x-independent Gaussian contribution, and 



SmfM - / d'-ep{\e\)SMF[M;fi~£] 



(6.37) 



where p(|e|) is the (assumed isotropic) single site distri- 
bution for the site disorder e, and Smf[M; — e] is the 
pure action (|6.25p with fj, ^ fi — e. It follows that the 
phase boundary is now defined by: 



d'^T / d'^ep{\e\)e 



-A(e,)C/o|r|^~ = --(/x-£)-^ 



(6.38) 

For Ci- = 1, it is assumed here that |m — ^1 < Uo/2 for 
all allowed values of e, otherwise the mod factor in (|6.33p 
enters. 

The result (I6.38[) . even for = 1, does not reproduce 
the result fJ.c{Jo', A) — Hc{Joi 0) — A (except at Jq — 0), 
where A is the maximum value of |e| supported by p, 
obtained from the rare region arguments in Sec. IIII Dl 
In fact, the integral (|6.38|) is smaller than it would be 
if all the weight of p were to lie at e = — A/i (whence 
\fj,\ + A), implying that Jq.c (ImI) is larger, leading 
to enlarged Mott lobes: Hc{Jo', A) > fic{Jo', 0) — A. 

In the presence of disorder, (|6.38[) implies that the 
Mott lobes shrink, but in a nontrivial way that depends 
on the detailed shape of p, not just its support. This is 
an artifact of mean field theory, which couples all sites 
equally, whereas the rare region argument result relied 
on local interactions in spatially separated, noninteract- 
ing droplets. Moreover, the coupling constant Jq/N in 
(|6.2ip is scaled by the total number of sites. Therefore, 
even if one were to view all sites with values of Si in some 
small region of size (Ae)'^^ as a droplet, its critical point 
would have to be scaled by the relative number of sites 
in the droplet, i.e., by l/p{\e\){AeY^ , in order to obtain 
a consistent result. This roughly explains the form of 
(I05D . 

For er = I, the Mott lobes still disappear entirely for 
A > C/o/2. However, for < 1, the fact that the Mott 
lobe extends to arbitrarily large means that it will sur- 
vive for arbitrarily large A as well, including unbounded 
disorder. 



A/„ = --m 



^g-/3'=-(27rn//3-iM) 

II 



2Uo 2/3^- 
1 



^In 



27r 



■In 



1 



0, er <l 



/2 



(6.40) 



in which, once again, ^o(m) — ~^ ^ (nonzero only 

for — 1) is the integer that minimizes {I — ^/Uq)'^, 
and the last line follows for /3 00. The derivative 
with respect to fi therefore recovers the Mott lobe integer 
density p = lo(fJ,) for Cr — 1 and p = for e,- < 1. The 
vanishing compressibility follows immediately. 

In the presence of disorder, loifJ- — s) is un- 
changed within a (smaller) Mott lobe, so AfQ^^{fi) = 
/ d'^^ep(|£|)A/o(/x — e) = A/o(/.t) is also unchanged, and 
identical results follow for the density and compressibil- 
ity. 



5. Phase diagram for finite range hopping 

In both pure and disordered cases, the fact that Mq « 
y^—TQ/uo is nonzero immediately outside of the Mott 
phase boundary, with unchanged mean field critical be- 
havior, implies that there is no glassy phase in mean field 
theory. Only for finite-ranged hopping will a finite width 
Bose glass phase appear between the Mott and superfluid 
phases, as in Fig.[TI^b). 

In considering the interpolation between small Et- and 
Et- = 1 it is clear that it makes sense only to consider the 
zero density Mott lobe. The renormalization group treat- 
ment in Sec. IVI Dl will therefore focus on the vicinity of 
the commensurate transition at = 0. The mean field 
calculations give one some confidence that the critical 
behavior near this point will vary continuously with er, 
and that extrapolations from small values to the physical 
value = 1 will be at least qualitatively valid. 



4- Mott phase compressibility 

Computation of the density and compressibility within 
the Mott lobe requires the M-independent part of the free 
energy /o in (|6.26p . Using (|6.25p and (|6.28p one obtains 
/o = foo + A/o, where 



/oo = -^In 

I r d'- 

" '2 J (2^ 



D(j){T)e 

27rC/o 



■In 



(6.39) 



C. Droplet model for general e-r 

In order to gain insight into the nature of the Bose 
glass phase for general er, we next generalize the droplet 
model of Sec. IIIII The identical rare region argument as 
for er — I implies that the Mott phase boundary occurs 
at /i±( Jo, A) — /i±( Jo, 0) =F A, where A is the maximum 
supported value of e^. For /i| > |^±(Jo, A)| there will 
. We consider then the effective action (|6.9I) in a finite 
superfluid droplet, with spatial volume V, in which is 
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assumed constant in space [compare p.30p ]: 



S. 



(1) 

eff 



V / d'^T 



ip ■ Vr(t 



..41) 



in which, as before, p, k are the bulk density and com- 
pressibiHty if the droplet were extended into an infinite 
homogeneous medium. The usual 27r-periodic boundary 
conditions on (j) apply. 



1. Droplet density and compressibility 

Let us first generalize the free energy calculation, 
p.32p . By performing a calculation identical to (|6.39p 
and (|6.40p . the generalization of (|3.36p in the vicinity of 
chemical potential ^to then reads 



(6.42) 



■In 



in which the functions /o and /oo are defined analogously 
to (jXM)) and (0351): 



foo{no) 



(2^ 



■In 



27r 



(6.43) 



is the free energy associated with the ^k"(Vt-0) term 
[compare (j6.39p ]. and 

/o(|/x|) = /o(lMol) - - ImoI) - - |/xo|)^ 

(6.44) 

is the analytic background free energy in the neighbor- 
hood of /io- In the limit /3 — > c», for e^. = 1, one recovers 
the discrete particle addition result ()3.39p and ()3.40p . For 
Er < 1 only the 1 = term in (|6.42p survives, and one 
obtains 



/ — /o + .foo + TTZo ■ 



(6.45) 



Noting that = k1, one obtains 



Pdrop(Mo) = - 



„ 1 (dp 
= p" p° = 0. 



(6.46) 



The superfluid droplet therefore has vanishing density, 
and hence, like the Mott phase, is incompressible, for 
arbitrary p and V. We conclude that for e,- < 1 both the 
Mott phase and the Bose glass phase are incompressible, 
and K (and p) will become nonzero only as one crosses 
into the superfluid phase. 



2. Droplet temporal correlations and superfluid 
susceptibility 

One must therefore seek a different measure to distin- 
guish the Mott and Bose glass phases for Et- < 1 . We ex- 
amine, therefore, the temporal correlation function and 
superfluid susceptibility. 

By following through the derivation (|3.4ip - (|3.43p for 
er < I, analogous to the survival of only the 1 = term 
in the free energy series (j6.42p . the (3 ^ oo limit for the 
temporal correlation function leads to the simple result 



,i[0(-r)-0(x')] 



= Go(x-x')e-'"(^-^')/'" (6.47) 



in which 



Go(x) 



' IkV 



(6.48) 



is the Gaussian correlation at p = 0, with coefficient 
A(tr) defined by ((6301)— compare ((02l) . 

Clearly the Go factor, which decays faster than ex- 
ponentially, dominates the asymptotic behavior at large 
|t — t'I, but for large kV this may not occur until |x — t'| 
is extremely large. To quantify this, we compute the 
droplet averaged correlation function 



-'^i^-r I dVdnp{V,K)Go{T) (6.49) 



where p{V, k) is the probability density for a droplet of 
size V and bulk compressibility k. The last line holds in 
the vicinity of the Mott lobe boundary, where p ~ kqAp, 
where Ap = /x — /Xc(Jo,A) is the deviation from the 
phase boundary, and ko(Jo) is the compressibility just 
inside the background superfluid phase — compare (j3.3ip . 
Therefore p-T / k w Ap ■ t is approximately independent 
of K. 

For large |t|, the integral (|6.49p is dominated by large 
volumes for which piV^n) ~ e^^/^°(''\ where Vo{k) is 
a characteristic scale for droplets with background com- 
pressibility n. A steepest decent evaluation of (|6.49p then 
leads to 



Gp{T) 



r/2 



(6.50) 



where To(Jo) = ['«oK)('«o)/4A(er)]^^'^~'^^', and various 
non-exponential prefactors have been dropped. The 
|t|/to)^~'^^/^ exponent reproduces (|3.20p at Ap = and 
Ct — 1. This is quite a remarkable result: Although the 
temporal correlations decay at infinity within any given 
droplet, the large droplet tail of the distribution causes 
the droplet averaged correlation function to diverge in 
directions such that Ap - t — > —oo. This strongly distin- 
guishes the Bose glass phase from the Mott phase. 
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A related quantity is the superfluid susceptibility, 
Xs = [^)^ ^ / ^'''/ (6-51) 



The saddle point estimate of the contribution from a sin- 
gle droplet is given by, 



Xd: 



,.op(F) « V / d^-TG<f\T) 



^(l-e^)Ko" 



(6.52) 



which grows faster than exponentially in V, and is there- 
fore dominated by large droplets. The total susceptibility 



(6.53) 



therefore diverges, providing a more direct signature of 
the Bose glass phase. 



3. Droplet model conclusions: er < 1 vs. er = 1 

The droplet model exhibits some strong differences in 
the physics of the Bose glass phase, all related to the 
physics of superfluid droplets, between e,- = 1 and e^- < 1 
that may raise questions about the smoothness of the 
limit e-r ^ ^■ 

The vanishing density and compressibility in the Bose 
glass phase implies that there should be no analytic back- 
ground contribution to k for e,- < 1. Thus (|6.17p . with 
some exponent oi(er) < 0, is expected to describe the 
critical behavior of the full compressibility, not just its 
singular part. This means that the analytic contribution 
to the density and compressibility must appear discon- 
tinuously at e,- = 1. 



Due to droplet incompressibility, the mechanism lead- 
ing to a divergent susceptibility for < 1 is also rather 
different than that for Et- = 1. As discussed in Sec. lIII CTSl 
for the latter, the temporal correlations decay as a slow 
1/t power law — see (|3.24p and (|3.48[) — due to the spec- 
trum of finite sized droplets that are very close to adding 
or giving up an extra particle. For < 1, the long-time 
correlations are dominated by large droplets, and it is 
the resulting stretched exponential behavior that gener- 
ates the divergent susceptibility. 

These strong differences do not require that critical 
exponents, such as a, be discontinuous, but one could 
certainly imagine some kind of nonanalytic behavior, as 
when one approaches an upper or lower critical dimen- 
sion of a transition. We shall see below that, at least, 
the qualitative features of the critical behavior do be- 
have as functions of Cr in the expected fashion. Since 
the expansion is around e,- = 0, it has absolutely nothing 
to say about possible singularities at e^- — 1. However, 
the results do lend some support to it as an analytic tool 
for exploring some features of the Bose glass-superfluid 
transition in higher dimensions where numerical results 
are not yet available. 

D. Renormalization group calculations 



1. Replicated Lagrangian 



In this final section, we revisit the renormalization 

1 — I 

group calculation carried out in Ref. 16, but now with 
explicit attention to issues of particle-hole symmetry. To 
this end, we use the standard replica trick on Cc in (EH) 
to average over the disorder,— and obtain the replicated 



Lagrangian, C 
ber of replicas, and 



where p — > is the num- 



OL — 1 



A3[^*V,V - golV-l'Kx, t) ■ [V/V,V - go|V'l'](x,r') 



(6.54) 



(6.55) 



r 



In the end we will take er — Cx — 1, but in setting 
up the RG calculation it is useful to leave them as free 
parameters. For simplicity we have taken (5g(x) and 
i5f(x) = 5r(x) — (5g(x)^ to be independent Gaussian ran- 



dom fields with 

[5g(x)],, = 0, [^g^(x)5g.(x')],, = A,<5(x-x')<5^. 

[<5r~(x)]^, = 0, [<5r~(x)<5r~(x')],,-A.<5(x-x'). (6.56) 

The analysis in Ref. [l^ was carried out with go finite 
and Ag = 0, but A^ > 0. However, considerations of 
particle-hole symmetry now lead us to seek fixed points 
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FIG. 9: Vertices corresponding to u, and Ag, where 
are replica indices, and the u vertex also enforces energy con- 
servation, LOl — LU2 -\- 1^3 ~ i^i ~ 0. 



with go = and > 0. 
In Fourier space we have 

p 

- It. [ I I I '^(kl - k, + k3 - k4) 

Q, ^—1 "^ki •^k2 ks k4 

X ■(/'a(ki,'^)'0a(k2,w)'(/;;^(k3,w')V'/3(k4,'^')> 

(6.57) 

where = A,. + gpAg. Nominally, by naive power 
counting, it would appear that the leading term at low 
frequencies should be the A^ term, and one might expect 
the other two frequency dependent terms to be strongly 
irrelevant. This will turn out to be true for small Cr, 
where naive power counting is almost valid. However, 
because u) and uj' refer to different replicas, a and /3, 
these terms break particle-hole symmetry, and are not 
as strongly irrelevant as one might expect (as compared 
to, for example, corrections to the coefficient, u), 
and, in fact, will be seen to become relevant beyond a 
critical value of e,-. 

We shall also begin by setting go — 0. By power count- 
ing, the go • V'*Vt-V' term (focused on in Ref.[2) again 
appears strongly relevant compared to the CtIVt-V'Pi a-^d 
the Aggo(?/;*Vr'!/')|'!/'P appears strongly relevant com- 
pared to the Ag(V'*VTV') • (V'*VtV') term. Again, al- 
though relevant at small e.^, there is a (different) critical 
value beyond which these terms become irrelevant. 



2. Recursion relations 

We shall perform a standard Wilson momentum shell 
renormalization group transformatior^ in which succes- 
sive shells in fc-space are integrated out. For each such k 
all frequencies, w, are integrated out. Since the frequency 




FIG. 10: Diagrams that contribute to the propagator renor- 
malization. 




FIG. 12; Diagrams that contribute to the renormalization of 
Ar. 
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FIG. 13: Diagrams that contribute to the renormahzation of 



is unbounded, the Brillouin zone is really a hypercylin- 
der. After each integration, we rescale k and u) in order 
to maintain the same wavevector cuttoff, k{^. The spin 
rescaling factor and the dynamical exponent, z, are de- 
termined in the usual way by setting and e,- equal to 
unity. 

In Fig. [9] are shown the four basic diagrammatic ver- 
tices corresponding to u, A,., and Ag, while in Figs. fTDHT^ 
are shown the lowest order diagrams that then contri- 
bution to the renormahzation of the propagator and of 
the vertices themselves. Note that there are a number 
of "missing" diagrams not included because they do not 
contribute in the replica limit, p — > 0. We obtain, then, 
in a straightforward way the recursion relations: 



dr 
di 

du 

lu 

dAr 
"df 

dA, 



2(m + l)u 2Ar ^, 9 

1 + r 1 + 7' " 



eu - 2(m + 4)u^ + 12uAr + 0{u^, . 



(6.58) 



(6.59) 



dl 



(e + er)Ar + 8Al - 4(m + l)uAr + 0{u^, . . .) 

(6.60) 

(6.61) 



^ - Ag(e + e^ + 10A,. -2Ag-2) 
: A,A,-2A^, 



where m is the number of boson species (m = 1 phys- 
ically), f = ro/kl, u = Kdu/'i, Ar = KdAr, Ag = 
k\KdAg are appropriately rescaled by the cutoff, and 
Kd = 2/(47r)''/2r(d/2) is (27r)-'* times the area of the 
unit sphere in d-dimensions. The conditions e^, = e,- = 1 
lead to the identifications 



Z=l + Ar 



0. 



(6.62) 



Dirty bosons fixed point 



Note that Ag does not enter any recursion relations 
except its own at this order. If one sets Ag = one ob- 
tains the usual Boyanovsky-Cardy lowest order recursion 



relationsfi^ with fixed point 

_^ 3me + {5m + 2)er 



A* = A* = 

9 ' ' 



, ^ e + 3er 
8(2m- 1) ' - 4(2m- 1) 
(2 - m)e + {■m + 4)e^ 



8(2m- 1) 



(6.63) 



correct to linear order in e and e^- For sufficiently small 
e, Cr we see that A* < and, consistent with the previ- 
ous power counting estimates, this fixed point is stable 
against the perturbation Ag. However, for 

13m + 16 

this fixed point becomes unstable. Setting to = 1 and 
e = (so that e,- = 1 corresponds to d = 3) we find 



29' 



(6.65) 



which is actually quite small, and is therefore a poten- 
tially reasonable estimate. Using (|6.6ip and (|6.62p at 
Ag = 0, one may write Ag = 2z — d + SA^. The last 
term shows the rather large deviation from the naive re- 
sult Ag = 2z — d, discussed in Sec. IVI A 31 The latter 
would have led to the estimate = |, which is uncom- 
fortably close to unity, considering how poorly controlled 
this expansion is for larger trr^ 

For Et- > e^, there is now a new stable finite Ag fixed 
point, 

a; = i(e + .. + 10A:-2) = i|!i±^(.,-.^), (6.66) 

which bifurcates continuously away from the random rod 
fixed point. The corresponding dynamical exponent, 

,-1 , A* I j.._ (^ + 4)^+(7^ + i0K ^ 

z - i + A,, + Ag - ^^2^ _ , > e„ 

(6.67) 

is substantially larger than the random rod value. For 
TO = 1, e = 0, and Cr = 1 one obtains z = which 
should be considered a very crude extrapolation. The 
"thermal" eigenvalue, determining v, is 



2-2(to + 1)u* + 2A;! 



= 2 



(to + 4)e + (7m + 10)6^ 
8(2to- 1) 



(6.68) 



while, as stated above, rj = 0. These results are both 
unchanged from their random rod values at this order. 



4- Relevance of particle-hole symmetry breaking 

Finally, let us include the go term. To linear order the 
flow equation for go — I go I is found to be 



dgo 
dl 



go[l + A,- Ag]. 



(6.69) 
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FIG. 14: Proposed behavior of the random rod (RR), com- 
mensurate (statistically particle-hole symmetric) dirty bo- 
son (C), and incommensurate (fully particle-hole asymmetric) 
dirty boson (IC) fixed points as functions of e-r. The Ar axis 
has been suppressed. Here G and Gb are commensurate and 
incommensurate Gaussian fixed points. The commensurate 
fixed point bifurcates away from the random rod fixed point 
at ~ At e-r = Cri — |, C and IC merge and at er = 1, 
C is the stable fixed point that describes the physical dirty 
boson problem. 



Consistent with the power counting estimates, go is 
strongly relevant for small tr, and is always relevant at 
the random rod fixed point, Ag = 0. However, go be- 
comes irrelevant at the dirty boson fixed point for 



4(2™ - 1) 



3(1 



(6.70) 



For m — 1 and e = one obtains e^-i = |. This is again 
an uncontrolled estimate, but does indeed indicate that 
the statistical particlc-holc symmetry is restored prior to 

tr - 1. 

For e-T < €ri there is a new incommensurate fixed point 
at nonzero qq. In order to locate it one should choose z 
to keep rather than e-r, fixed during the RG fiow. The 
flow equation for A„ then becomes 



dl 



= Ag[e + e^-4 + 8A^]+8A2, 



(6.71) 



while the remaining flow equations are identical to those 
in Ref. [l^. The flxed point found in that work did not 
account for Ag, but we see that if A* = 0{e,er) then 
A* = 0(e^, ecr, e?), so the results given there are indeed 
correct to 0(e, Cr). However, as er grows, so does A*, 
eventually leading to the merging with the dirty boson 
flxed point. 

The technical problems encountered in Ref. iS are also 
overcome in our analysis. Specifically, the |Vt-V'P term, 
which was also previously ignored, has the flow equation 



der 



(6.72) 



implying that it is of the same order as A* at the fixed 
point. A nonzero Ct fixes the convergence problems, and 
allows one to remove an unphysical frequency cutoff wa- 
Finally, to see how the two fixed points merge, we write 
down flow equations in the intermediate region, e^i > 
tr > f^, where one must consider both Cr and = k^go- 
We choose z so that Ct -f go = 1 remains fixed. The flow 
equation for go is then 



dm 
dl 



(2 + 2A^ - z)go 



with 



z = (l-l-A^Ag)-f (l + A^)go. 
At the fixed point we therefore find 



1 + a: 



9o 



1 



(6.73) 



(6.74) 



(6.75) 



which vanishes precisely when go becomes irrelevant at 
the dirty boson fixed point. 

The entire proposed fixed point structure as a func- 
tion of Cr is summarized in Fig. I14p . For small er the 
unstable random fixed point and stable filly particle-hole 
asymmetric fixed points exist. For e^ < e,- < Et-i there 
are three fixed points, with the new commensurate dirty 
boson fixed point being more stable than the random rod 
fixed point, but less stable than the asymmetric fixed 
point. Finally, for e-r > Cri the incommensurate fixed 
point merges with the dirty boson fixed point, which is 
then completely stable. This provides a detailed scenario 
by which statistical particle-hole symmetry is restored . 
We caution, however, that due both to the uncontrolled 
nature of the double e-expansion at the dirty boson fixed 
point, and the special nature of e^ = 1, extrapolation 
of these results to er = 1 should be treated as, at best, 
qualitative estimates. The general scenario we propose, 
however, seems very natural and illuminating. 



APPENDIX A: FUNCTIONAL INTEGRALS 
BASED ON THE TROTTER DECOMPOSITION 

Given a Hamiltonian, 7i, thermodynamics is obtained 
from the partition function Z = tr [e~^''~^] . The Trotter 
decomposition involves identifying a convenient complete 
set of states {|a)}, writing e~^'^ = [e"^"^^] , where 
M = P/At, and inserting the states |q;) between each 
element of the product: 



EE--- E("oi 



|ai)(ai|e" 



"2 



X ... X (aA/-i|e 



(Al) 



Often one can decompose H = Hq+Hi, and choose the 
\a)'s to be eigenstates of TCq. This is especially convenient 
when H = TiUqi}, {pi}] is written in terms of a set of N 
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canonically conjugate positions, ^i, and momenta, pi, and 
takes the special form 

n^Ho[m]+ni[{p.,}] (A2) 

in which Tii is a quadratic polynomial in the pf. 

Hi = i ^ S^JP^PJ + ^ Up,. (A3) 

The |a)'s are then taken to be product eigenstates of the one obtains 



positions qi: \q) 



Ljqi) with qi\q) 



qi\q)- We also define product eigenstates of the mo- 



menta, \p) = \pi, . 



,Pn) 



LiIPj), with p,\p) = p,\p). 



From the canonical commutation relations, [qi,Pi 
which in turn imply 

[e'^'^\Pj] = -\e'^^^6,_ 
we immediately infer that 



e-^'P^q,) 



\qj + X) 
\p,+X). 



(A4) 



(A5) 



From this we obtain the wavefunctions 



{q^\p^) = {q^ - 0\e"''^'\p^) = 6^?'^ ' (g, = 0|e^f'«'|K = 0) 
= ^e^'-'^ (A6) 

where = ili = ^\Pi = 0) normalizes the wavefunction. 
Note that if the qi are periodic variables, defined only 
modulo 2tt, say, then \qi) = \qi + 2tt) are identified. This 
requires the same periodicity of the wavefunctions, and 
determines the allowed eigenvalues, pi — which therefore 
must be integers in this case, and one obtains Af = \f2Ti 
as well. 

With the shorthands q = {tji}, p — {pi}, etc., we may 
now compute, for small At, 



^^|g-ArWo(9)g-ArWi(p)|^/^ 



-ArHoiq) 



g-ArWi(p)|^/^ 



-A 



-^«(l)J2e-^^'''^P\q\p){p\q') 



-ArWo(g) 



(A7) 



and we are now left only with computing the inverse 
Fourier transform of the Gaussian function e^^'^^^'^-'. 
To do this, we first complete the square: 



ui = 'Y{S^^)ijtj] u ^^SijVj. (AS) 



We specialize now to the case of integer pi. Using the 
formula 



^ = / dp, Y 5{p, " rii) 

-/ — OO „ _ 



rii — — oo 



/oo 
dpr 
-oo ^ . — 



(A9) 



— — OO 



.»E.Pi(9>-<z:)p-AT-Hi(p) 



= J2 [ -^e*^^P'(9--«'+2-'"-)e-^-«i(p) (AlO) 

■m 



dp 



g-l^r J SijPiPj+iY.i Pi(qi-q'i+2-n 



AA(At) ^ 

X e"5S7 E,.3(S"').i(9i-g;+2™.)('3j-9j+2Tmi) 

where, in the second equality, we have changed variables 
to Pi =pi + Vi, and A/'(Ar) = ^det(ArS'). 

Now consider the limit Ar — > 0. For given {qi] only 
a single term in the m-sum will survive, namely that 
which minimizes the exjDonent. Furthermore, only if 
qi — q'i + 27rmi = 0{At^) will the term contribute to 
the path integral, IjAip . Therefore, modulo 27r, qi{T) be- 
comes a continuous function in the limit At — > 0, and 
Qi ~ Qi + ^nnii — > — At^^. Clearly we will nearly al- 
ways have TTii = 0, with = ±1 every so often when 
the periodic boundary conditions are enforced (for exam- 
ple, when qi = 2n~ and g- = 0^, or vice versa). When 
At — > we may equivalently define a continuous func- 
tion qi^r), taking arbitrary real values, and then sum 
over all boundary conditions qiir) — qi{0) + 2'Kmi, with 
mi running over all integers. Thus, we finally have, as 
At 0: 



AA(At) 



-iArE,.,(S- 



3 9i«J 



(All) 



X e 
1 



Af{AT) 



,-iATE,.j(S"'),j('ji+iti)(<j,+it,) 



and 



Vqij) exp 



dr 



Ho[g(T)] 

(A12) 
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where J 'Dq{T) is a functional integral over all paths with 
a uniform (Wiener) measure. 

Equation (jA12p now leads directly to the path integral 
representation for the Josephson Hamiltonian (j2.3p . with 
q replaced by </> and p replaced by n. The term H{q) is 
just the Josephson cosine coupling term. 

If, instead of canonical coordinates, Ti, = Ti,[a\a] is 
written instead in terms of raising and lowering opera- 
tors, another convenient complete set of states is that of 
the coherent states, 



{a\a') 



E 

oo 

E 

n=0 



rt), a\a) — a\a) 



(A13) 



Thus, for any normally ordered operator, 0{a\a), we 
have 



{a\0{a\a)\a') = 0{a* , a'){a\a') 
and hence for At — > 



{a\e 



^ArW[at,a]i /\ 



(a|(l- ArH[a^a] 



(A14) 



(A15) 



Recognizing that 



E 



At 



(A16) 



where the periodic boundary conditions in imaginary 
time have been used, we arrive at (|2.ip (with tp replacing 
a) when At — > 0. 



APPENDIX B: DUALITY TRANSFORMATIONS 

In order to obtain a model amenable to analysis in one 
dimension one must derive a dual representation for the 
Josephson Lagrangian, (|2.3p . In Rcf. 6 this was done 
in a somewhat ad hoc fashion using the Haldane repre- 
sentation for one-dimensional bosons. Here we perform 
the duality transformation directly, on a variant of the 
Josephson Lagrangian (j2.3p . in a much more transparent 
manner, following closely the analogous derivation for the 
two-dimensional AF-model and its associated Kostcrlitz- 
Thouless transition^ 

The transformation is best carried out in discrete time. 
The continuous time limit is mathematically well defined 
but, as we shall see, physically less transparent: one runs 
into logarithmically divergent coupling constants, just as 
in the Trotter decomposition of the quantum Ising model 



in a transverse field.^ These are consequence of the usual 
exponential weighting times in continuous time, discrete 
state Markov processes. In order to avoid introducing 
the probabilistic formalism necessary for dealing with the 
continuum limit we shall maintain a discrete time vari- 
able. 

We begin by introducing the Villain, or periodic Gaus- 
sian form, of the AF-coupling:'^'^ 



^-/fo(l-cos(0)) 



E 



e 



Vq(4>,k) 



(Bl) 

which allows the duality transformation to be carried 
out exactly. We will consider only the case where in 
(|2.3|) is nearest neighbor, but possibly random. In (|B1|) 
Kq = JijAr for the bond (ij), and therefore already in- 
cludes the effect of discretizing the T-integral in (|2.3p . 
There are two limits in which K and Kq may be quanti- 
tatively compared. For large Kq the variable will have 
only small fluctuations about zero, and only the curva- 
ture, ATqi near the minimum of the cosine potential at 
= is important. In this limit only the m = term 
contributes to the Villain form, and the two potentials 
therefore match when K « Kq. Conversely, when Kq is 
small, fluctuates strongly, and many m contribute to 
the Villain sum. It is then convenient to use the Fourier 
representation^^ 



,Vo{(t>,K) _ 



/— — oo 



where now, for small K, only the I 
important. This yields 



(B2) 



0, ±1 terms are 



VQ{d^,K) 



1 



ln(27rii') -(-2e-^cos((/)), ^ 0, (B3) 



and hence the correspondence K w l/21n(2/i4ro) in this 
limit. Now, the continuum limit. At — s- 0, corresponds to 
Kq — > 0, and therefore (|B3p is appropriate. This yields 
K K, — 1/2 ln(Jy At) as At — the logarithmic be- 
havior alluded to above. We shall keep At small but fi- 
nite, setting aside the question of its optimal value, which 
would have to be addressed for quantitative comparison 
of the Josephson and Villain forms of our model. For our 
purposes it is important only that the discrete and con- 
tinuous time versions lie in the same universality class. 

Using (|B1[) we then define the Villain form of the 
Josephson Lagrangian, (|2.3p : 

m 

£j ^ ^ Knda.4>r-2TTm'if 



X (9,0,, -iz/, -27rmO,), (B4) 
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where da(t>r = 



"r+a 



r and nir = (to°, mj, . . . , mf) 
is a (d + l)-diinensional integer vector field defined at 
each space-time lattice site r = (i, t), and the index a = 
0, 1, . . . , d runs over the neighboring sites in the a = Xq 
direction, with xq — 'T'- Until further notice there is 
no restriction on the dimensionality, d. The parameters 
iff, Vij and Vi are, respectively, the Villain analogues of 
JijAr, UijAr and fiiAr in (|2.3p . 

Since e^-' is separately periodic in all the differences 
9a4>r we may write it as a Fourier series. 



where 
1 c 

so that 



2 ^ iff 



(B6) 



and the normalization is 



M= /det(27rV) J] 

V i,r,Q#0 * 



(B7) 

(B8) 



g — 7711;^— m /2K 



To derive (IB7[) we have used the identity 

^^-irnu^~^K{u+ivf _ ]_ 

J 27r V2^' 

(B9) 

and its higher-dimensional generalizations. The partition 
function now becomes 



(BIO) 



where V • n is the discrete space-time divergence, 

V.nr = 5]«-<_J. (Bll) 

a 

This formulation is entirely real, and is therefore a con- 
venient basis for Monte Carlo simulations of the dirty 
boson problem:" 

Let us now restrict attention to d = 1. One may then 
solve the constraint V • n = by introducing a dual lattice 
integer field, 5'r, such that 

rir - (V X S)r = (5R_i - 5r, Sr ~ 5R-f ), (B12) 

where the dual lattice bond connecting R — t to R is 
the one that cuts the real lattice bond connecting r to 
r -t- X, while that connecting R — x to R cuts the one 
connecting r to r -|- f, i.e., R = (i, T) = r -I- i(x -I- r). 
Thus the a-component of the discrete curl of a scalar is 



the difference between its values on the two dual sites 
that border the bond from r to r -I- a. The field S'r is 
defined uniquely up to an overall additive constant, and 
the constrained trace over the Mr is precisely equivalent 
to the free trace over the 5r. One therefore obtains 



_^ £.7[VxS] 



(B13) 



with 



A7[Vx5] = -i^^(5.5R)2 

R ^ 



I,J,T 



(B14) 



R 



Here daSn = 5'r — S'r^q, Ki is the Villain coupling 
on the real lattice bond that cuts (R — x, R), and simi- 
larly for i^i and Vij. Note that in this representation the 
Lagrangian is purely real and has a very natural clas- 
sical interpretation, namely that of a three-dimensional 
interface model. The field S'r represents the height of a 
surface over a two-dimensional plane. The first two terms 
in Cj determine the energy cost for steps in the r and 
X directions, respectively. In this case the energy associ- 
ated with steps in the r direction is random, but only in 
the spatial index. The last term represents a random tilt- 
ing potential which favors steps in the x direction with 
the same sign as vj. It is precisely this breaking of the 
symmetry between left and right steps that reflects the 
broken particle-hole symmetry in the original quantum 
Hamiltonian. Note that in this dual model the symmetry 
being broken is associated with parity (x —x) rather 
than time reversal (r — r). 

The more common sine-Gordon representation is ob- 
tained from (jBSp by softening the integer constraint on 
the 5r. Thus E7^=-oo - ir^dSnEZ=-.oHSn - 
/ir) is replaced by J dS-R,q{S-R), where q(t) is periodic 
with period one, and is peaked around t ^ 0. The sine- 
Gordon model (|5.ip results from the choice 



^2yo cos(27rt) 



(B15) 



where yo is called the fugacity. The integer constraint is 
recovered in the limit yo ~* oo- For small j/o (see below) 
this term may be obtained directly by including a term 



ln(yo 



R 



(B16) 



in the original Lagrangian, (jB4[) . The discrete curl of a 
vector field is a scalar field on the dual lattice obtained 
by summing the vector field around the dual lattice pla- 
quette. 



(V X m)R 



(BIT) 
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and is precisely the vorticity on that plaquette. 

Finally, the Coulomb gas representation is obtained 
either from the sine-Gordon representation by expand- 
ing the exponential in (|B15[) and integrating out the 
S'r, or from the discrete version (jB14p . by substituting 
J:Z=-oo e'''^'-^- for EZ=-oo S{Sn " hn): 



Z = tr'tr'^ 



tr' 



(B18) 



where [we include the term (jB15l) for completeness], 



R,R' 



(B19) 



R 



and 5rr' is the inverse of the quadratic form: 

(.G-')-R-R' = -^{dTdT'Sn-R') + ididrVH')STT'. (B20) 
Kj 



For diagonal Vij = VqSij and fixed Kj = Kq, Grr' 
is, modulo a trivial rescaling, the inverse of the two- 
dimensional lattice Laplacian, and yields the usual log- 
arithmic Coulomb interaction at large distances. So 
long as Vij is short ranged and Kj = Kq + SKj with 
SKj/Kj ^ 1, ^RR' will remain Coulomb- like at large 
distances. Note that Cc is once again complex, with dxi^ 
playing the role of complex offset charges. 



The sine-Gordon form yields the same Coulomb gas 
form (jB19[) . except that the values of Zr are restricted 
to 0, ±1 only. When yo is small, large values of /r are 
suppressed anyway, and the difference between (jBlSP and 
(|B16P is negligible. 



^ See, e.g., N. Read and S. Sachdev, Phys. Rev. Lett. 62, 
1694 (1989); A. V. Chubukov and S. Sachdev, Phys. Rev. 
Lett. 71, 169 (1993); J. Miller and D. A. Huse, Phys. Rev. 
Lett. 70, 3147 (1993). 

^ See, e.g., J. A. Hertz, Physica Scripta TIO, 1 (1985) and 
references therein. 

^ H. P. Wei, D. C. Tsui, M. A. Paalanen and A. M. M. 
Pruisken, Phys. Rev. Lett. 61, 1294 (1988); A. M. M. 
Pruisken, Phys. Rev. Lett. 61, 1297 (1988); S. C. Zhang, 
T. H. Hansson and S. Kivelson, Phys. Rev. Lett. 62, 82 
(1989); N. Read, Phys. Rev. Lett. 62, 86 (1989); X.-G. 
Wen and Y.-S. Wu, Phys. Rev. Lett. 70, 1501 (1993). 

* See, e.g., P. A. Lee and T. V. Ramakrishnan, Rev. Mod. 
Phys. 57, 287 (1985). 

^ See, e.g., (a) D. B. Haviland, Y. Liu and A. M. Goldman, 
Phys. Rev. Lett. 62, 2180 (1989); (b) A. F. Hebard and 
M. A. Paalanen, Phys. Rev. Lett. 65, 927 (1990); (c) M. 
P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990). See also A. 
M. Goldman and N. Markovic, Physics Today, November 
1998, pp. 39-44, and references therein. 

^ M. P. A. Fisher, P. B. Weichman, G. Grinstein and D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 

F. Pazmandi, G. T. Zimanyi and R. Scalettar, Phys. Rev. 
Lett. 75, 1356 (1995); K. Sheshadri, H. R. Krishnamurthy, 
R. Pandit, and T. V. Ramakrishnan, Phys. Rev. Lett. 75, 
4075 (1995). The mean field approximation fails to pro- 
duce a Bose glass phase except as a line separating the 
Mott insulating and superfluid phases. It does, however, 
produce a region of quasi-localized superfluid phase, with 
exponentially small couplings between well separated su- 
perfluid droplets, which should be the seed of the Bose- 
glass phase, and thus cross over to the truly localized non- 
superfluid phase, when fluctuations about the mean field 
ground state are taken into account. 

* J. K. Freericks and H. Monien, Europhys. Lett. 24, 545 



(1994); Phys. Rev. B 53, 2691 (1996); M. Niemeyer, J. K. 
Freericks, and H. Monien, Phys. Rev. B 60, 2357 (1999). 
^ Y. Tu and P. B. Weichman, Phys. Rev. Lett. 73, 6 (1994); 
Y. B. Kim and X. G. Wen, Phys. Rev. B 49, 4043 (1994). 
M. B. Hastings, Phys. Rev. B 64, 024517 (2001). 
^° L. Zhang and M. Ma, Phys. Rev. B 45, 4855 (1992); K. G. 
Singh and D. S. Rokhsar, Phys. Rev. B 46, 3002 (1992). 
Many authors have explored various models in one and 
two dimensions using quantum Monte Carlo techniques: 
See, e.g., (a) R. T. Scalettar, G. G. Batrouni and G. T. 
Zimanyi, Phys. Rev. Lett. 66, 3144 (1991); (b) W. Krauth, 
N. Trivedi and D. M. Ceperley, Phys. Rev. Lett. 67, 2307 
(1991). (c) M. WaUin, E. S. S0rensen, S. M. Girvin, and 

A. P. Young, Phys. Rev. B 49, 12115 (1994). (d) S. Zhang, 
N. Kawashima, J. Carlson and J. E. Gubernatis, Phys. 
Rev. Lett. 74, 1500 (1995); (e) F. Alet and E. S. S0rensen, 
Phys. Rev. E 67, 015701 (R) (2003); (f) N. Prokof'ev and 

B. Svistunov, Phys. Rev. Lett. 92, 015703 (2004). 

M. Makivic, N. Trivedi and S. UUah, Phys. Rev. Lett. 71, 
2307 (1993). There is, however, some argument over the 
interpretation of the unusual values of the exponents ob- 
tained in this work: see P. B. Weichman, Phys. Rev. Lett. 
74, 1038 (1995); N. Trivedi and M. Makivic, Phys. Rev. 
Lett. 74, 1039 (1995). 

J.-W. Lee, M.-C. Cha, and D. Kim, Phys. Rev. Lett. 87, 
247006 (2001); J.-W. Lee and M.-C. Cha Phys. Rev. B 72, 
212515 (2005). 

A. Priyadarshee, S. Chandrasekharan, J.-W. Lee, and H. 
U. Baranger, Phys. Rev. Lett. 97, 115703 (2006). 
(a) S. N. Dorogovtsev, Phys. Lett. 76A, 169 (1980); (b) D. 
Boyanovsky and J. L. Cardy, Phys. Rev. B 26, 154 (1982); 
Erratum, Phys. Rev. B 27, 6971 (1983); (c) I. D. Lawrie 
and V. V. Prudvikov, J. Phys. C 17, 1655 (1984). 
P. B. Weichman and K. Kim, Phys. Rev. B 40, 813 (1989). 
^"^ R. Mukhopadhyay and P. B. Weichman, Phys. Rev. Lett. 



42 



76, 2977 (1996). 

E. Altman, Y. Kafri, A. Polkovnikov, and G. Refael, Phys. 
Rev. Lett. 93, 150402 (2004); larXiv:0711.2070k l [cond- 
mat.dis-nn]. 

I. F. Herbut, Phys. Rev. B 57, 1303 (1998). 
1. F. Herbut, Phys. Rev. B 58, 971 (1998). The 1 + e expan- 
sion developed here consists of computing the RG flows in 
1-D, and inserting the engineering dimensions to general- 
ize the recursion relations to higher dimensions. However, 
in the absence of a form for the Hamiltonian for noninte- 
ger d > 1, there is presently no rigorous support for this 
approach. 

Further neighbor interactions substantially increase the 
complexity of the phase diagram in the absence of the 
random site energies, ei. One can, in principle, generate 
Mott insulating phases with arbitrary rational densities 
("charge density wave" states). The superfluid transitions 
from these states are surprisingly complex: for a review, 
see L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and 
K. Sengupta, Prog. Theor. Phys. Suppl. 160, 314 (2005). 
Generically, only the integer fillings are stable against small 
disorder since the fractional fillings necessarily break the 
lattice translation symmetry, leading to multiply degener- 
ate ground states related by a discrete translation. It is 
not hard to see that arbitrarily small random ei will al- 
ways generate rare regions where it is energetically favor- 
able to form two such states with a domain wall between. If 
one allows further neighbor hopping matrix elements, Jij, 
with various signs, one can also generate supersolid phases 
which break both lattice translational symmetry and XY- 
phase symmetry, i.e., superfluid charge density waves: for 
some recent work see, e.g., P. Sengupta, L. P. Pryadko, 

F. Alet, M. Troyer, and G. Schmid, Phys. Rev. Lett. 94, 
207202 (2005); G. G. Batrouni, F. Hbert, and R. T. Scalet- 
tar; and references therein. 

See, e.g., A. L. Fetter and J. D. Walecka, Quantum Theory 
of Many Particle Systems (McGraw HiU, 1971), Chaps. 10, 
14. For a more modern view, see P. B. Weichman, Phys. 
Rev. B 38, 8739 (1988). 

See, e.g., J. V. Jose, L. P. Kadanoff, S. Kirkpatrick and D. 
R. Nelson, Phys. Rev. B 16, 1217 (1977), and references 
therein. 

This idea has been used to explain the vanishing of the 
Hall conductivity at magnetic fleld-tuned superconduct- 
ing transitions: See M. P. A. Fisher, Physica A 177, 553 
(1991); A. T. Dorsey and M. P. A. Fisher, Phys. Rev. Lett. 
68, 694 (1992). 

In d = 1 arbitrarily weak hopping disorder in the spin- 
i XXZ chain with vanishing axial magnetic fleld^'^ leads 
to an insulating phase consisting of bound (ferromagnetic) 
singlet pairs — a random singlet glass: see D. S. Fisher, 
Phys. Rev. B 50, 3799 (1994). In the ID boson prob- 
lem this means that the superfluid transition takes place 
at finite Jo,c(/ifc) > even at half-integer fiUing;'^^ and 
Fig. [He) must be adjusted accordingly, with a line of ran- 
dom singlet glass (identified by both an infinite superfiuid 
susceptibility,"^ and an infinite compressibility, k — > oc^^) 
a,t jl — jlk on the interval < Jo < Jo,c{P-k)- 
In the limit of small Jo/Uo, only site occupancies of k 
and -|- 1 contribute between Mott lobes, and (|l.ip or 
(|1.2|) may be mapped onto a spin-i XXZ model Ti.! — 

^I^(i,j> JiA^i^j in which h oc fi-{k + 
^)Uo vanishes at half-filling [see, Ref.-^, as well as M. E. 



Fisher, Rep. Prog. Phys. 30, 615 (1967)]. If /i = 0, the ener- 
getics, for d > 1^- favors spins aligned in the plane. So long 
as the model is not singular, in the sense that there is fi- 
nite probability po that a given Jij vanishes, the model will 
have long range superfluid order, -i/jo = [{d'f + icrf )]av 7^ 0. 
Moreover, in this Jo/Uo limit, Jo is the only energy 
scale, and the quantum state must be independent of Jo, 
showing that superfluidity survives for arbitrarily small Jo . 
If Po lies above the bond percolation threshhold, the lat- 
tice will break up into finite, non-communicating droplets, 
and bulk superfluidity is suppressed entirely. For nonzero 
h, roughly speaking, those sites where h/ Ji lies above some 
threshold of order unity, where Ji — ^ j J^j 1 align 
with h along z (or —z, for h < 0). For sufficiently large \h\, 
depending on the precise distribution of Jij, the z-aligned 
spin clusters will percolate, and superfluidity will be de- 
stroyed. The resulting state is the spin-^ analogue of the 
Bose glass phase, and the remaining isolated planar or- 
dered clusters (those with anomalously small h/ Ji) are the 
superfluid droplets of Sec. lIIII The two Mott phases corre- 
spond to \h\ sufficiently large that all spins are a-aligned 
(or anti-aligned) . 

For a detailed discussion of the onset transition, see, e.g., 

P. B. Weichman, Phys. Rev. B 38, 8739 (1988). 

D. S. Fisher and M. P. A. Fisher, Phys. Rev. Lett. 61, 1847 

(1988). 

^® P. B. Weichman and R. Mukhopadhyay, Phys. Rev. Lett. 
98, 245701 (2007). 

There are, however, also not so obvious fractional time 
dimension generalizations of (|2.H) : see Ref . [l6i . 
The superfluid susceptibility, defined by Xs ~ 
{d^o/dhs)hs=o, where hs is a symmetry breaking field 
entering via coupling —hs "^^i&i + a|) added to p.ip . or 
— hs cos(0i) added to 1)1. 2[) . is related to the full space- 
time pair correlation function via Xa = J d'^x J dTG{:x.,T). 
Although G(x, r) decays exponentially in space, with 
a finite correlation length, its more singular temporal 
behavior leads to distinct behaviors in the different 
insulating phases. The stretched exponential decay (|3.20|) 
is integrable, so Xs is finite in the random rod glass (as 
well as in the Mott phases). However, the power law decay 
p.24p in time leads to a logarithmic divergence of Xa, so 
that the Bose glass phase is characterized by vanishing 
superfiuid order, but an infinite superfluid susceptibility.— 
^2 R. B. Grifliths, Phys. Rev. Lett. 23, 17 (1969). 

M. E. Fisher and V. Privman, Phys. Rev. B 32, 447 (1985). 
The random site energies may be thought of as locally re- 
ducing the local chemical potential fi-Si = iJ, + A — {A+ei) 
from a uniform value fi + A. A free energy convexity argu- 
ment shows that that such a reduction can only reduce, or 
at best maintain, the particle density. Since, by assump- 
tion, fi + A < /i+(Jo) lies within the Mott phase, the den- 
sity will remain fixed at the Mott integer value. 
D. S. Fisher and P. C. Hohenberg, Phys. Rev. B 37, 4936 
(1988). 

^® There are dynamical scaling breakdown issues at the finite 
T lambda transition [described by the classical Model F 
equations: P. C. Hohenberg and B. I. Halperin, Rev. Mod. 
Phys. 49, 435 (1977)] as well, where the finite k argument 
now yields z — d/2. However, violations are possible (and 
believed to occur in d = 3) where there exist two different 
dynamical exponents zi, zi satisfying 21+22 = d. Only the 
mean {z\ -\- Z2)/2 = d/2 enters the corresponding hydro- 



43 



dynamic correlation function [V. Dohm, Phys. Rev. B 44 
2697 (1991)]. 

There are actually two distinct sets of scaling functions 
and amplitudes in l|4.7|l . one for 5 > and one for S < 0. 
For notational simplicity we shall not make this distinction 
explicit since ps (and k in the random rod problem case) 
are nonzero only in the superfluid phase. 
No nonuniversal amplitude A is required in this formula- 
tion due to quantum hyperuniversality: see K. Kim and P. 
B. Weichman, Phys. Rev. B 43, 13583 (1991). 
The definition of the superfluid density (|4.8|l in terms of 
finite size scaling of twisted boundary conditions, was first 
proposed by M. E. Fisher, M. N. Barber and D. Jasnow, 
Phys. Rev. A 8, 1111 (1973). Amazingly, the last section 
of this paper, containing a discussion of how finite size 
scaling might lead to violations of the Josephson relation, 
proposes a mechanism very similar to ours: anomalously 
strong boundary condition dependence of te, the distance 
from the critical point, on the twist wavevector, fco. If tg — 
t ^ ko they find v — 1 — a. Here we find an even more 
anomalous linear dependence of Sg — 5 on ujo, leading to 
Vt — —a. We are unaware of any classical model where the 
original proposed violation occurs, but a candidate would 
presumably be a transition between two different phases 
with finite superfiuid density. 

A famous theorem [J. T. Chayes, L. Chayes, D. S. Fisher 
and T. Spencer, Phys. Rev. Lett. 57, 2999 (1986); Com- 
mun. Math. Phys. 120 501 (1989)] requires that Ufs > 2/d, 
where ffs is a correlation exponent defined through finite 
size scaling, and it is commonly assumed that under most 
conditions that = f^. However, it has been argued more 
recently [F. Pazmandi, R. Scalettar, and G. T. Zimanyi, 
Phys. Rev. Lett. 79, 5130 (1997)] that Vfs generally places 
no constraints on the value of u. Nevertheless, analytic re- 
sults in d = 1 (Ref. !6l) and quantum Monte Carlo data in 



d = 2 (Ref. [Tib are consistent with this inequality. 
Irrelevance of go means that it drops out completely suffi- 
ciently close to the dirty boson critical point, ruling out in 
yet a different way, any scaling (|4.7|l of shifts go — > go — i^^o 
with ^T- 

T. Giamarchi and H. J. Schulz, Europhys. Lett. 3, 1287 
(1987); Phys. Rev. B 37, 325 (1988). 

In this interface picture, the singlet glass phase at half 
filling2^ corresponds to sufficiently small Ko = [-ft'/]av in 
which the hierarchy of different values of Kj (organized 
from smallest to largest) generates an unusual hierarchy 
of long-range spatial correlations where all sites divide up 
into (possibly well separated) pairs. Although the individ- 
ual tilts diSf^j^T), diS(ji^T) at two such sites /,/' fluctuate 
strongly, as a function of T, between and 1, the sum 
diSi^j^T) + diS(ii^T) remains close to unity — the path in- 
tegral analogue of the quantum state -^(|0, 1) -I- |1, 0)).— 
This pairing reduces the interface tilt fluctuations at large 
distances, reducing the roughness below that of the super- 
fluid phase. Since our main concern in this paper is with 
the neighborhood of integer flUing, we shall not pursue the 
physics of this singlet glass any further here. 
I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Se- 
ries and Products, (Academic Press, New York, 1980). 
See, e.g., A. Aharony, Phys. Rev. B 12, 1038 (1975); 
see also in "Phase Transitions and Critical Phenomena," 
edited by C. Domb and M. S. Green (Academic, London, 
1976), Vol. 6, Chap. 6. 

See, e.g., K. G. Wilson and J. Kogut, Physics Reports 
12C, 75 (1974). 

See, e.g., M. Suzuki, Prog. Theor. Phys. 56, 1454 (1976). 
See, e.g., E. S. S0rensen, M. Wallin, S. M. Girvin and A. 
P. Young, Phys. Rev. Lett. 69, 828 (1992). 



